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We calculate the self-force experienced by a point scalar charge q, a point electric charge e, and a 
point mass m moving in a weakly curved spacetime characterized by a time-independent Newtonian 
potential $. We assume that the matter distribution responsible for this potential is bounded, so that 
$ ~ —M/r at large distances r from the matter, whose total mass is M; otherwise, the Newtonian 
potential is left unspecified. (We use units in which G = c = 1.) The self- forces are calculated by first 
computing the retarded Green's functions for scalar, electromagnetic, and (linearized) gravitational 
fields in the weakly curved spacetime, and then evaluating an integral over the particle's past world 
line. The self-force typically contains both a conservative and a nonconservative (radiation-reaction) 
part. For the scalar charge, the conservative part of the self-force is equal to 2£q 2 Mr /r 3 , where £ is 
a dimensionless constant measuring the coupling of the scalar field to the spacetime curvature, and 
f is a unit vector pointing in the radial direction. For the electric charge, the conservative part of 
the self-force is e 2 Mf/r 3 . For the massive particle, the conservative force vanishes. For the scalar 
charge, the radiation-reaction force is ^q 2 dg/dt, where g — — is the Newtonian gravitational 
field. For the electric charge, the radiation-reaction force is |e 2 dg/dt. For the massive particle, the 
radiation-reaction force is —^-m 2 dg/dt. Our result for the gravitational self-force is disturbing: a 
radiation-reaction force should not appear in the equations of motion at this level of approximation, 
and it should certainly not give rise to radiation antidamping. In the last section of the paper we 
prove that while a massive particle in a vacuum spacetime is subjected only to its self-force, it is 
also subjected to a matter-mediated force when it moves in a spacetime that contains matter; this 
force originates from the changes in the matter distribution that are induced by the presence of the 
particle. We show that the matter-mediated force contains a radiation-damping term that precisely 
cancels out the antidamping contribution from the gravitational self-force. When both forces are 
combined, the equations of motion are conservative, and they agree with the appropriate limit of 
the standard post-Newtonian equations of motion. 



I. INTRODUCTION AND SUMMARY 



Motion of massive bodies in strong gravitational fields 

The problem of determining the motion of N bodies subjected to their mutual gravitational interactions has been 
a center piece of general relativity from its inception. Work on this started shortly after the introduction of the 
theory, and in 1938, a firm formulation of the equations of motion was given by Einstein, Infeld, and Hoffmann [Q, 
who provided post-Newtonian corrections to the Newtonian equations of motion. Higher-order corrections were later 
added and work continues today, with Damour, Jaranowski & Schafer 00), de Andrade, Blanchet & Faye 

and Pati & Will |n| currently computing corrections of third post-Newtonian order. A technical review of this field 
of research, as of 1987, can be found in Ref. JO]; the early history was described by Havas |l2| . 

While the equations of motion for slowly moving bodies in a weak gravitational field are now very well understood, 
the same cannot be said of the fast motion of a massive body in a strong gravitational field. The prototypical problem 
here is to determine the motion of a structureless particle of mass m in the gravitational field of a much more massive 
black hole. While it is known that the motion is geodesic in the limit m — * 0, the issue is to determine the corrections 
to geodesic motion that appear when m in nonvanishing. An important effect that must be introduced is the loss of 
orbital energy and angular momentum to the gravitational radiation emitted by the moving particle; this is effected 
by an appropriate radiation-reaction force. Another important effect is a conservative correction to the equations of 
motion, which disappears when m — > 0. Schematically, therefore, the equations of motion will take the form 

= /self, (LI) 

where u a is the particle's four-velocity in the background of the massive black hole, u a .^uP its acceleration (the 
semicolon denotes covariant differentiation in the black-hole metric), and / s " lf is the self- force. This contains both a 
conservative and a dissipative (radiation-reaction) component, and it scales as m 2 in the small-mass limit. (Through- 
out the paper we work in relativistic units, with G and c both set equal to unity.) 
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A useful way to look at Eq. (1.1) is to imagine that the motion is actually geodesic in a spacetime that is not 
the background spacetime of the massive black hole, but the perturbed spacetime that contains the particle as well. 
The geodesic motion in the perturbed spacetime can be expressed as a forced motion in the background spacetime, 
and this gives rise to Eq. fll.lf ). The implementation of this idea, however, presents some difficulty: If the particle is 
pointlike, the perturbation diverges at the position of the particle, and the geodesic equation is not defined on the 
world line. The perturbation must then be decomposed into a part that is singular but does not influence the motion 
of the particle, and a smooth remainder that is entirely responsible for the self- force Jl3| . This decomposition requires 
great care. 

In simple situations, the orbital evolution of a particle emitting gravitational waves can be determined without 
the involvement of a self-force. If the black hole is nonrotating, the orbital evolution is determined by energy and 
angular-momentum balance |14 lq |: The rate at which the radiation carries energy and angular momentum both to 
infinity and across the event horizon matches the rate at which orbital energy and angular momentum is lost by the 
particle. This information alone is sufficient to solve for the motion, at least when the time scale for orbital evolution 
is long compared with the orbital period. If the black hole is rotating, the statement continues to be true provided 
that the orbit is either equatorial or circular []l6|-^l| . When dealing with generic orbits around a rotating black hole, 
however, the loss of orbital energ y a nd angular momentum no longer constitutes sufficient information — the rate of 
change of the "Carter constant" [0 must also be determined, and this requires the involvement of a self- force |^0| . 



Sources of gravitational waves for LISA 



There is a practical necessity for the computation of the gravitational self-force. Solar-mass compact objects in 
highly eccentric, nonequatorial orbits around rapidly rotating, massive black holes (in the range between 10 3 and 
10 8 solar masses) emit gravitational waves that will be measured by eventual spaceborne interferometric detectors 
p3[ , which operate in the low- frequen cy b and (in the range between 10~ 4 and 1 Hz). One such detector, the Laser 
Interferometer Space Antenna (LISA) [Q , has been selected as one of three "Cornerstone missions" in the "Horizon 
2000+" program of the European Space Agency. A possible involvement by the National Aeronautics and Space 
Administration — now appearing likely after the publication of the Taylor-McKee decadal survey (2f| , in which LISA 
is listed as a high-priority, moderate-size mission — would increase the likelihood that this detector will be deployed 
in the not-too-distant future. A realistic projection places the launch in the neighborhood of 2010. 

The detailed modeling of these sources of gravitational waves for LISA, to the extent that templates could be 
provi ded for data analysis, will require a detailed computation of the orbital evolution. This, in turn, requires that 
Eq. (fTl]) be evaluated and solved for generic orbits around a Kerr black hole. For this we need a practical way of 
computing the gravitational self-force. 



Gravitational self-force 

This issue was taken up recently by Mino, Sasaki, and Tanaka po| , who calculated the self- force acting on a point 
particle of mass m moving in an arbitrary background spacetime with metric g a p] they assume that the metric satisfies 
the Einstein field equations in vacuum. Their calculation is based on a careful analysis of the perturbation field in 
the vicinity of the particle, the (gravitational) perturbation being created by the particle itself. It involves a careful 
decomposition of the perturbation into a part that is singular at the particle's location but does not affect its motion, 
and a part that is well-behaved and governs the motion. Their expression for the self-force was then reproduced by 
Quinn and Wald p7fl , on the basis of a different approach involving a comparison between self-forces acting in two 
different spacetimes. Both teams found the following equations of motion: 
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P - /cxt - T m (^ + fL + /grav (1-2) 



Here, u a is the particle's four-velocity in the background spacetime, f" xt is an external force acting on the particle, 
and an overdot indicates differentiation with respect to proper time r. The gravitational self-force is given by 

/ g Q rav = -2m 2 f (2G V„' ;7 - G w / + u a G^, 5 u s ) u *Wt/ dr'. (1.3) 

The four-index object G a ^ , s , (x, x') appearing inside the integral is the retarded gravitational Green's function [p8[, 



whose precise definition is given by Eq. (3.25) below (our normalization differs from the conventional choice by a 
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factor of 4); the unprimed indices refer to the field point x (which is identified in the integral with the current position 
of the particle), while the primed indices refer to the source point x' (identified with the particle's past position). The 
integral extends over the past world line of the particle, from r' = — oo to (almost) the current time, r' = r~ = r — e, 
where e is very small and positive P^j. The integration is cut short to avoid the singular behavior of the Green's 
function as x 1 approaches X] it involves only the smooth part of the Green's function, which is often referred to as its 
"tail part" |p9fl . In a situation in which the particle is subjected only to its self- force, the equations of motion reduce 
to mu a .puP — fg Tav - Because the self- force scales like m 2 , we recover geodesic motion in the limit m — > 0. 

The difficulty of evaluating Eq. (O) resides with the computation of the gravitational Green's function, which 
would be very hard to carry out for an arbitrary spacetime. There is hope for progress, however, if the spacetime 
possesses symmetries, such as staticity and spherical symmetry in the case of a Schwarzschild black hole, or stationarity 
and axial symmetry in the case of a Kerr black hole. In such cases, a plausible method of computation would be 
based on a separation-of-variables approach, and the modes of the Green's function could be obtained fairly easily 
pfj| . But this approach is not free of difficulties: While the individual modes of the Green's function stay finite 
and continuous as x approaches x' (though they are discontinuous in their first derivatives), the sum over modes 
does not converge. Essentially, this is because the modes do not distinguish between the singular and smooth parts 
of the Green's function; they contribute to both, and the singular behavior of the Green's function gives rise to a 
divergent sum. Recently, Ori, Burko, and Barack (3l]-|37]] have devised a way to regularize the mode sums, so as 
to extract from them a meaningful expression for the self-force. They have applied their technique to a number of 
simple situations involving scalar, electromagnetic, and gravitational radiation. Similar regularization methods were 
proposed by Lousto as well as Nakano, Mino, and Sasaki |3!^|l(|. It appears likely that in the near future, this 
method will be used successfully to solve more complicated problems, including the prototypical problem mentioned 
previously. 



Electromagnetic self-force 



The derivation of Eqs. (1.2) and (1.2) by Mino et al. |26J], and Quinn & Wald ]27|] , is based on methods previously 
introduced by DeWitt and Brehme [jllfl. These authors calculated the self- force acting on an electrically charged 
particle moving in an arbitrary spacetime with metric g a p. As corrected by Hobbs }42^| , the equations of motion of a 
charged particle are 



= /« t + ~ + u a up)fi t + l -e 2 (R a p uP + u a R p ^u^ + (1.4) 



where e is the particle's electric charge, f" xt an external force acting on the particle, R a p the spacetime's Ricci tensor, 
and 

L = -e 2 f (G a Y0 - G^; a )u^' dr' (1.5) 

J — oo 

is the electromagnetic self-force. The two-index object G a p,(x,x') is the retarded electromagnetic Green's function 
pl[ , whose precise definition is given by Eq. (3.16) below. In the absence of an external force, and in a region 
of spacetime empty of matter, the equations of motion reduce to mu a .pvP = f" m . In flat spacetime, on the other 
hand, /™ n = because the Green's function contains only a singular part; the smooth, or tail, part vanishes. In flat 
spacetime, therefore, we recover the Abrahams-Lorentz-Dirac equation 



V 3 = /« t + ~ (6 a p + u a up)fi t . (1.6) 

The practical computation of the electromagnetic self-force presents the same technical difficulties as in the gravita- 
tional case. 



Scalar self-force 

The self- force acting on a particle with scalar charge q was recently calculated by Quinn [jl6). In this case, the 
equations of motion are 
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X - /Set + + u a up)fL + ^q 2 (R a p u? + u a R^u f3 u-<) + /« calar , (1.7) 



3m v p M/JCXI 6 

where the self-force is now given by 



Ciar = 1 2 {a aP + u a uP) / G, dr', (1.8) 



in which the retarded Green's function G{x,x') is a scalar; this is defined by Eq. (3.1C) below. In the absence of an 



external force, the equations of motion reduce to mu a pU /3 = / s " a i ar , and this is the simplest realization of the class 
of equations that govern the orbital evolution of a particle moving in a strong gravitational field. The computation 
of the scalar, electromagnetic, or gravitational self-force involves solving for a retarded Green's function, and then 
performing an integration over the particle's past world line. Because there is only one component to the s cala r 



Green's function (as opposed to 16 for the electromagnetic function, or 100 for the gravitational case), Eq. fll.q ) 
captures the essence of the problem and avoids many technical complications. For this reason, it has been the starting 
point of many recent investigations ^ |37],|3^,[l(],[l7j . 

Weak-field limit 

The computation of the various Green's functions is a very difficult undertaking even when the spacetime possesses 
many symmetries. This is, however, a very tractable task when the spacetime is only weakly curved. In such a 
situation, the metric itself can be calculated perturbatively as an expansion about flat spacetime, and the same 
technique can be employed to find the Green's functions and compute the self-forces. Our purpose with this paper is 
to do just that. We will consider spacetimes for which the metric can be expressed as 

ds 2 = -(1 + 2$) dt 2 + (1 - 2<&){dx 2 + dy 2 + dz 2 ), (1.9) 

in which &(x,y,z) <C 1 is a generic Newtonian potential satisfying Poisson's equation, V 2< 1> = 47rp, where p is the 
mass density. We will keep $ unspecified throughout the paper, but assume that it is small everywhere. We shall 
also assume that the mass distribution is bounded, so that <E> behaves as —M/r at large distances r from the center 
of mass; M — J pd 3 x is the total mass, and r 2 = x 2 + y 2 + z 2 . Various aspects of those spacetimes are discussed in 
Sec. II. In Sec. Ill and IV we calculate the scalar, electromagnetic, and gravitational Green's functions to first order 
in <f>. In Sec. V we use the Green's functions to calculate the scalar, electromagnetic, and gravitational self-forces. 

Our motivation for carrying out such a (long) computation is to provide a useful check on the formalism, especially 
in its gravitational formulation. Admittedly, the weak-field limit is not by itself very interesting, but understanding 
this limit may well be a necessary step toward understanding the strong-field behavior of the gravitational self-force. 
So while we cannot hope to learn anything new from such a weak-field computation, we can use our knowledge of the 
post-Newtonian equations of motion 1 1 1 to produce a nontrivial check of the self- force method; the results found 



here should of course match the results from the literature. As we shall see, we have been successful at producing this 
check: the method works. 

DeWitt and DeWitt 

We are not the first to perform such a check on the formalism. In 1964, DeWitt and DeWitt Gq] calculated the 



electromagnetic self- force, as given by Eq. fll.q ), for the weakly curved spacetime of Eq. (L9). They considered only 
the special case $(cc) = —M/r, where r = \x\, but noted that their results generalize to arbitrary potentials by 
superposition. (We use a bold-faced symbol to denote a three-dimensional vector living in flat space.) They found 
that to leading order in a weak-field, slow-motion approximation, the spatial components of the self-force are 



o M . 2 9 dq 

f em (x)=e 2 —r + -e 2 -£, 1.10 

r } 3 dt 



where f = x/r and g(x) = — V$ is the Newtonian gravitational field. In this limit, the equations of motion 
■mu a pvP = f" m become 

d 2 z 

m-j^ =mg(z) + f cm (z), (1.11) 
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where z(t) is the trajectory of the charged particle. The first term on the right-hand side of Eq. ( L.1C ) is a repulsive 
correction to the local force of gravity; it agrees with the weak-field limit of a result by Smith and Will [|49|, who 
calculated the force required to hold a charged particle in place in the (exact) field of a Schwarzschild black hole. The 
second term is the usual expression for the radiation-reaction force experienced by a charged particle subjected to an 
external force mg E3M4I]. 

The result of Eq. ( [LIO; ) is remarkable because this entirely local expression for the self-force derives from a nonlocal 
formulation involving the entire past history of the charged particle; in the weak-field limit, any trace of nonlocality 
is lost. The result is remarkable also because while the expression for the radiation-reaction component of the 
self- force is the expected ^e 2 d g/dt . it is very hard to see how such a simple result could ever follow from such a 
complicated expression as Eq. (1.5). We notice that this result for the radiation-reaction force can be obtained from 
two very different approaches. In a flat-spacetime point of view, the equations of motion for a charged particle are 
mu a pvP — fc Xt + §( e2 / TO )(^ Q /3 + u a u p)fe Kt ; in the slow-motion limit, the external force is mg, and we recover the usual 
result. In a curved-spacetime point of view, on the other hand, there is no external force and the equations of motion 
are mu a l3 u 13 = f" m ; in the weak-field, slow-motion limit, we also recover the usual result. This agreement is necessary 
on physical grounds, but it could ha rdly h ave been anticipated on the basis of a quick inspection of Eq. ( [1.5] ); the 
calculations required to produce Eq. ( 1.1C| ) are very involved. 

In their paper pSfl , DeWitt and DeWitt gave the following physical picture for the self- force of Eq. ( l.lCj ). This 
picture relies heavily on their expression for the Newtonian potential, $ = —M/r, which corresponds to a point mass 
M located at the origin of the coordinate system. 

Because the fundamental expression for the self-force, Eq. (1.5), involves only the smooth, or tail, part of the 
electromagnetic Green's function, there is a priori no local contribution to this force coming from the world-line 
integral. Indeed, DeWitt and DeWitt found that the quantity inside this integral vanishes in the immediate past of 
the event (t,x), which corresponds to the particle's current position. In fact, it vanishes until the time delay, t — t' , 
becomes equal to r + r' , the distance between the central mass M and the particle's current position x, added to 
the distance to the particle's old position x' . At t — t' — r + r' , the electromagnetic Green's function undergoes an 
abrupt change of behavior, and the self-force is entirely due to this sudden transition. (We expand our discussion of 
this point in Appendix B.) The physical effect responsible for the force, in DeWitt and DeWitt's view is a signal 
that originates from the particle at an earlier time t', propagates toward the central mass M at the speed of light, 
bounces off the central mass, and comes back to the particle at the current time t. Although the self-force is nonlocal, 
Eq. (1.1C) involves the conditions at the current time only. This is because the time delay in not noticeable at the 
level of approximation maintained in the calculation. To leading order in a weak-field, slow-motion approximation, 
the electromagnetic self-force appears to be entirely local. 

The physical picture described in the preceding paragraph appeared to us to be slightly suspicious. The sudden 
change in the behavior of the electromagnetic Green's function at t — t' = r + r' is dictated by the presence of the 
point mass M at x = 0, which mediates the interaction between the charged particle and its earlier self. The region 
of spacetime near the central mass is therefore seen to have an important effect on the Green's function, and it is this 
effect — the sudden change — which apparently gives rise to the self-force. But $ = —M/r cannot be considered to 
be small in this region of spacetime, and the perturbative method of calculation of the Green's function must come 
into question. While there is little reason to doubt the validity of Eq. (1.10), it appeared to us that its derivation by 
DeWitt and DeWitt Gq] left room for improvement. 



This paper 



Part of this paper is concerned with providing a derivation of Eq. ( 1.10 ) that is free of this criticism. Throughout the 
paper we assume that &(x) is everywhere much smaller than unity, and reduces to —M/r only far away from the mass 
distribution. Building upon DeWitt and DeWitt's work H], we introduce (in Sees. Ill and IV) techniques that allow 
us to compute the electromagnetic Green's function for such generic New tonia n potentials. In Sec. V we recalculate 



the electromagnetic self- force for these potentials, and we reproduce Eq. ( 1.1C ). While our result is compatible with 
the physical picture suggested by DeWitt and DeWitt, our derivation shows very clearly that the result is largely 
independent of the conditions near the mass distribution. We consider this improvement on the original derivation to 
be a significant contribution to this field. Our derivation is also more transparent, in the sense that the computational 
labor involved is much reduced compared with the original calculations of DeWitt and DeWitt. 

The techniques introduced in Sees. Ill and IV allow us to calculate all three types of Green's functions — scalar, 
electromagnetic, and gravitational. Those are used in Sec. V to calculate the respective self-forces. Part of the paper 
is therefore devoted to the calculation of the scalar and gravitational self-forces in the weakly curved spacetimes of 



Eq. (1.9). To the best of our knowledge, such calculations have not yet been presented in the literature. 
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Scalar self-force 



The computations are simplest for the case of a scalar charge, and we obtain 



,M 



T scala 



1 2 dg 



(1.12) 



where £ is a constant measuring the coupling of the scalar field to the spacetime curvatu re; a precise definition is given 
by Eq. ( p/7| ) below. The equations of motion for the scalar charge q are identical to Eq. (1.11), but with the scalar self- 
force replacing the electromagnetic self- force. For minimal coupling (£ = 0), the self- force is entirely dissipative, and 
we recover the expected result, ^q 2 dg/dt; this is the radiation-reaction force experienced by a scalar charge subjected 
to an external force rag. For minimal coupling, the self-force acting on a stationary charge vanishes; this agrees with 
Wiseman's result |47| for the force required to hold a scalar charge in place in the exact field of a Schwarzschild black 
hole. If £ > 0, the conservative part of the self-force is repulsive. If £ < 0, it is attractive. This result for the scalar 
self- force is very similar to Eq. ( 1.1 Op , and we notice that the highly nonlocal expression of Eq. (1.8) has managed to 



produce the expected local res ult. Once more the dual point of view holds: We can adopt a flat-spacetime point of 
view, set / s " a i ar to zero in Eq. (1.7), and get the correct r esul t by equating the external force to mg; or we can adopt 
a curved-spacetime point of view, set f^ xt to zero in Eq. (1.7), and get the correct result by evaluating the world-line 
integral of Eq. (|J). 



Gravitational self-force and matter-mediated force 



The computations required for the gravitational self-force are the most involved, but here we also face a serious 
technical problem. We have stressed the importance of working with a Newtonian potential that is everywhere small. 
In fact, the perturbative calculation of the Green's function requires that the deviations of the metric g a p with 
respect to the Minkowski values r\ a $ be everywhere small. To accommodate this requirement, it is necessary that 
the spacetime contain matter: all vacuum solutions to the Einstein field equations representing an isolated massive 
object necessarily possess regions of strong curvature. Such strong-field regions would affect the Green's function 
in a way that cannot be predicted by a perturbative calculation; the entire method would fail, even if the Green's 
function were to be evaluated only in the weak-field region of the spacetime. Thus, the global character of the Green's 
function, together with the weak-field limitations of our calculational methods, dictate that we work with spacetimes 
that contain matter. 

Th e p roblem resides with the fact that the equations of motion for a massive particle, as given by Eqs. ([L 



and (1.3), are restricted to vacuum spacetimes. The Mino-Sasaki-Tanaka-Quinn-Wald equations of motion [|26|j27| 
are therefore not directly suited to a weak-field calculation. An extension to spacetimes containing matter must be 
produced. We do this in Sec. VI, proceeding in two steps. 

First, we incorporate the modifications to the gravitational Green's function that come from the presence of Ricci- 
tensor terms in Green's equation. These modifications affect the gravitational self-force, but they do not alter Eq. ( |l.3| ) 
if the particle is restricted to move in a vacuum region of the spacetime. Second, we consider how the presence of 
matter modifies the equations of motion. Because the background stress-energy tensor T af3 depends on the metric, it 
necessarily suffers a perturbation when the massive particle is inserted in the spacetime. Physically, this corresponds 
to the fact that while the object of mass M — the star — is at rest when the particle is absent, it must move in the 
presence of the particle: both objects move around a fixed center of mass |)(|. This perturbative motion of the star 
produces an additional metric perturbation, over and above the perturbation directly associated with the particle's 
own stress-energy tensor. This must be accounted for, and we shall see that its effect is to modify the right-hand side 



of Eq. (1.2) by a term which we call the "matter-mediated" force. While the self- force can be thought of as a 
direct action of the particle on itself, this additional force can be thought of as an indirect action mediated by the 
presence of matter: the particle acts on the star, and the star acts back on the particle. 

The gravitational self-force is computed in Sec. V, and this calculation incorporates the effects of the matter on the 
retarded Green's function. We obtain 

/grav = __ m _ (1.13) 

and according to the naive equations of motion md 2 z/dt 2 = mg + / grav , the gravitational self- force does work on the 
particle at an average rate dW/dt — +^-m 2 \g\ 2 ; it therefore gives rise to radiation antidamping [pl[ . Notice that the 
self- force represents a correction of 1.5PN (post-Newtonian) order to the Newtonian equations of motion ]5^ . Notice 
also that the dual point of view seems at work also here: our expression for the "radiation-reaction force" could be 
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calculated on the basis of Eq. ([[j) by adopting a "flat-spacetime point of view" in which /" av would be set to zero 
and the "external force" equated to mg. 

The remarkable conspiracy that makes the tail integral reproduce the external-force term in the equations of motion 
is seen to be at play in all three cases. While the agreement between the flat-spacetime and curved-spacetime points 
of view is quite necessary in the case of the scalar and electromagnetic self-forces, it is decidedly disturbing in the 
gravitational case. How can we understand this result? 

The answer obviously comes from the matter-mediated force, which we define and compute in Sec. VI. For this we 
obtain 



/mm = m5g + lPN 
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,dg 

dt ' 



(1.14) 



where the first term represents the change in the star's Newtonian gravitational field associated with its motion around 
the fixed center of mass, the second term is a post-Newtonian correction to the Newtoni an fo rce mg, and the third 
term is a radiation-damping term that precisely cancels out the antidamping force of Eq. ( 1.13| ). (Such a cancellation 
was noticed a long time ago by Carmeli ]53|] , in the context of a very different f ormu lation of the equations of motion.) 

A precise expression for the matter- mediated force can be found in Eq. ( 6.54 ) below. When it is substituted, 
together with the self-force, into the equations of motion, we find that they take the form 



fz_ 

dt 2 



M 



5to 



Av-V)p- p-(Av~ 3V)v + 4(p • v)V 



0(2pn,to 2 ,M 2 ), 



(1.15) 



which contains no trace of a radiation-reaction force. Here, z(t) gives the position of the particle, and v(t) = dz/dt 
is its velocity; the star moves on a trajectory Z(t) with a velocity V(t) = dZ /dt. We use p = z — Z to designate the 
separation between the two objects, and p = \z — Z\. The equations of motion for the star are 



d 2 Z 

dt 2 



m 

7^3 



l + 2v 2 



3 I p- v 



p ~ 3(p • v)v \ + 0(2pn, to 2 , mM). 



(1.16) 



Our calculations cannot reproduce all the terms in the Einstein- Infeld-Hoffmann equations of motion jljjsij: Eq. ( 1.1 5| ) 
omits terms that are quadratic in m and M, while Eq. ( 1.16| ) neglects terms of order m 2 and mM. The formalism we 
use in this paper is based on linear perturbation theory, and it would be incapable of producing corrections of order 
to 2 . On the other hand, the corrections of order mM and M 2 are within the reach of the formalism, but in order 
to produce them we would need to go beyond the weak-field approximation considered in this paper. Within these 
limitations, however, we have complete agreement between the calculations presented in this paper and the standard 
post-Newtonian treatment of the two-body problem. 



Organization 



The technical part of the paper begins in Sec. II with a detailed discussion of the metric of Eq. (1.9), and a 
derivation of the scalar, electromagnetic, and gravitational wave equations for the weakly curved spacetimes described 
by this metric. In Sec. Ill we introduce the two-point functions A(x,x') and B(x,x'), and show how the scalar, 
electromagnetic, and gravitational Green's functions can be obtained from them by acting with differential operators. 
In Sec. IV we introduce methods to compute the two-point functions, and evaluate them in interesting limiting cases. 
[In Appendix A we expand our discussion of A(x,x'), and in Appendix B we provide an explicit computation of the 
two-point functions for the special case <& = — M/r] In Sec. V we use our preceding results for the Green's functions 
to compute the scala r, elect romag netic , and gravitational self- forces for the spacetimes of Eq. (1.9); our results were 
quoted in Eqs. (1.12), (1.10), and (1.13), respectively. Finally, in Sec. VI we introduce and compute the gravitational 
matter-mediated force, and show that it cancels out the antidamping force calculated in Sec. V. 



II. WEAKLY CURVED SPACETIME 



A. The spacetime and its geometric quantities 

The spacetimes considered in this paper have a metric given by 



7 



ds 2 = -(1 + 2$) dt 2 + (1 - 2<$>){dx 2 + dy 2 + dz 2 ), 



(2.1) 



in which <!>(a:) is a Newtonian potential, a function of the spatial coordinates x satisfying Poisson's equation, 

V 2 $ = 47rp, (2.2) 



where p(x) is the mass density. By virtue of Eq. (2.2) and the fact that 1 everywhere in the spacetime, the 

metric of Eq. (2J) satisfies the linearized Einstein field equations. Throughout the paper we will work consistently to 
first order in $. If the metric were representing the gravitational field of a point mass M located at the origin of the 
coordinate system, then p = MS(x) and $ = —M/r, where r — \x\. We will not need to adopt this particular form 
for the Newtonian potential (whi ch w ould violate the condition that $ be small everywhere) , and we shall work with 
generic potentials satisfying Eq. ( |2.2|) . We will, however, assume that 3> becomes equal to —M/r, with M = J pd 3 x, 
far from the matter distribution. Notice that we do not allow $ to depend on time. 

Introducing t a = dx a /dt as the timelike Killing vector of the spacetime, we write the metric tensor as 

9a!i = Va0 ~ 2$Xa/3, Xq/3 = ?7a/3 + 2t Q i/3, (2.3) 

where 7] a p = diag(— 1, 1, 1, 1) is the Minkowski metric, which was used to lower the index on t a . Throughout this 
section we will continue to lower and raise indices with the Minkowski metric, unless it is otherwise indicated. Notice 



that Xa/3 = diag(l, 1, 1, 1). From Eq. ( j2~3[ ) it is easy to derive ^f^g =1-2$ and g a/3 = r] al3 + 2$x Q/3 - These results 
hold to first order in $. 



The metric of Eq. ( |2.3| ) comes with the following Christoffel symbols: 

K P = -X^.P - X^,a + Xa^- (2.4) 

As a consequence of the relation <&. a t a = 0, they satisfy ^"^r^ = 0. The Riemann tensor is given by 

<*/3 7 = ~X^, a p + X^, ay ~ + X al ^p- (2.5) 

Contracting over the first and third indices gives the Ricci tensor: 

R a p = XotpV®, (2.6) 

where □$ = rj al3 <fr^ a p. Because $ does not depend on time, □<& = V 2 $ = 8 ab Q, a bi where x a denotes the spatial 
coordinates. The Ricci scalar is R = 2n$, and the Einstein tensor is G a ^ = 2t a t /3 n$. Equation ( ^.2j ) follows from 
the Einstein field equations with T a @ = pt a w; this represents a fluid of mass density p and negligible pressure at rest 
in the spacetime. 

B. Wave equations 

The field equation for a massless scalar field -0 in a curved spacetime with metric g a p is 

g aP t, af 3 - tR1> - -4tt M , (2.7) 

where the semicolon designates covariant differentiation compatible with the metric tensor, £ is an arbitrary constant 
measuring the coupling to curvature, and p is a given source term. It will be convenient to work with a densitized 
form for the field equation, 

£?[V] = -47r > /=fl/i, (2.8) 



where E[ip] stands for the left-hand side of Eq. ( J2 . 7| ) multiplied by \J—g. For the weakly curved spacetime of Eq. (p.3j), 
we find that 

E[ip] = OtP + ^t a t p i) >a p - 2£(D$)0, (2.9) 

where □ = r\ a ^d a dp is the flat-spacetime wave operator and, as was pointed out previously, = V 2< & = Airp. 

An electromagnetic field is represented by a vector potential A a which is here assumed to satisfy the Lorentz-gauge 
condition, A a . a = 0. Its field equations are 

gTA° - R%A» = -4-Kj a , (2.10) 
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where j a is a given current density. The wave equation in densitized form is 

B a [A]=-Airy/=gj a , (2.11) 

where stands for the left-hand side of Eg. ( p.lOj ) multiplied by \J—g. A straightforward computation reveals 

that for the weakly curved spacetime of Eq. ( |2.3| ) , 

E a [A] = UA a + 4SH?fA a tlu , - 2{ X a ^, v + - X^*'")*^ ~ [x°"*,^ + - X ^K- (2.12) 

Here, □ acting on the vector field A a still stands for the scalar wave operator of flat spacetime. 

A gravitational perturbation on a background metric g a p is described by a tensor 7 aj g, and the full metric is 
9a.fi = 5q/3 + 7a/3- By linearizing the Einstein field equations about the background, one obtains a wave equation 
for the perturbation. This equation takes a simpler form if it is expressed in terms of a "trace-reversed" field 
7 Qj a = j a p — \{g iWrs ( ll v)ga,pi which is then chosen to satisfy a Lorentz-gauge condition: j a/3 ,g = 0. (Here and below, 
indices on the perturbation field are raised and lowered with the background metric, and a semicolon designates 
covariant differentiation on the background.) With these choices, the Einstein field equations take the form of a wave 
equation for 7 Q/3 |2qj : 

g^r%u + iRf/r" + s«/r u = -ltesr*, (2.13) 

where 5T a P is the perturbation in the stress-energy tensor (so that the full stress-energy tensor is T Q/3 = T a/3 + ST a P, 
with the first term denoting the background's stress-energy tensor), R^avp is the Riemann tensor of the background 
spacetime, and S^ av p is constructed from the background's Ricci and Einstein tensors, 

S m vP — 2>R,i,(u9&)v — Rnv9a(3 — R9n{ a 9p)u — 2g^G a fj- (2-14) 



The parentheses around tensor indices indicate symmetrization with respect to these indices. Once Eq. (2.13) has 
been solved for the trace-reversed perturbation, the metric is recovered from the relation 7 Q( 3 = 7a/3 — \gaf3l, where 

7 = 5 Q/3 7a/3- 

The densitized wave equation is 

£ Q/5 [7] = -Vfay/^gST"?, (2.15) 



where E a P[y] stands for the left-hand side of Eq. (2.13) multiplied by y/—g. A rather long but straightforward 
computation reveals that in a weakly curved spacetime with metric (2.3), 

- U*(if#xv a , f f u + 2 7 Q/3 + UTlPnuvf") ■ (2-16) 
In the first term on the right-hand side, □ = ij^d^d^ is the scalar wave operator. 



III. GREEN'S FUNCTIONS 



A. Generating two-point functions 



The task before us in this section is the computation of retarded Green's functions for the wave equations considered 
in Sec. II — Eqs. ( |2.8[ ), (2.11), and ( 2. 15j ) . As we shall see, these Green's functions can be constructed by acting with 
differential operators on two generating two-point functions, A(x,x') and B(x,x'), defined by 



A(x,x') 



1 

2^ 



G ttat (x, x")<S>(x")G a&t {x'\ x') d 4 x" 



(3.1) 



and 



B(x,x')= I G tiat {x,x") P (x")G f{: * t (x",x')d A x" 



(3.2) 



Here, $ is the Newtonian potential of Eq. (2.1), p is the mass density of Eq. ((O), and G flat (ir, x' ) is the retarded 
Green's function of the flat-spacetime wave operator, 
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(3.3) 



this satisfies nG flat (x, x') = — AnS^x — x'), where 8±(x—x') is a four-dimensional ^-function, equal to S(t— t')5^{x— x'). 
This method to calculate the Green's function originated in the work of DeWitt and DeWitt Q|, who used it to 
compute the scalar and electromagnetic Green's functions for the special case $ = —M/r. Here we generalize their 
method to arbitrary potentials, an d to t he g ravitational Green's function. 

Inside the integrals of Eqs. (3.1) and ( |3.2| ), the function G flat (ir, x") picks out the past light cone of the point x, 
while G flat (x", x') picks out the future light cone of the point x' . (These light cones are those of the fiat background 
spacetime.) The intersection of the two light cones defines a closed two-surface S, and we see that A{x, x') is the 
average of the Newtonian potential over S, while B(x, x') is the average of the mass density. Both two-point functions 
are zero if x and x' are spacelike related, because the surface of intersection then disappears. 

Averages of derivatives of the Newtonian potential can be obtained by differentiating A(x, x 1 ). For example, 



and 



1 

2^ 



-!- / G &at {x,x")[d a ^(x")]G &a \x" ,x') d 4 x" = {d a +d a ,)A{x,x') 



G &a \x,x")[d a „^{x")]G Rat {x",x')d 4 x" = (daff + dcp + dcfl, +d a ,f 3 ,)A(x,x'). 



(3.4) 



(3.5) 



We use the notation d a /f — df /dx a , d a pi f = d 2 f /dx a dx 13 , etc., in which the tensor index (either a, a', or a") 
indicates wit h res pect to which variable (between x, x' , and x") the function / is being differentiated. 

Equations (3.4) and ( [3~5| ) are easy to establish. (In the following we keep the notation simple and drop the factor of 
1/2tt in front of the integrals). First we write d a A = J (d a G)QG dx" , and use the fact that G flat (a:,a;") depends only 
on the difference between x and x" . This allows us to express the right-hand side as — J (d a "G)&G dx" . Integrating by 
parts, we obtain + J G(d a »<S>)G dx" + J GQ(d a >'G) dx" . In the second term, the derivative op erat or can be switched 
to — dgi, which can be taken o utsi de the integral. The second term becomes — d a >A, and Eq. (|3.4D follows. To derive 



Eq. (|J) we differentiate Eq. (fj) with respect to x^ , and go through the same procedure 



Because the Newtonian potential does not depend on time, Eq. (3.4) implies 



d t A(x,x') + d t >A(x,x') = 0. (3.6) 

Anot her useful identity results from taking the trace of Eq. ( |3.5| ). Using Poisson's equation □$ = V 2 $ = 4irp and 
Eq. (3.2), we obtain 



-(□ + D')A{x, x') + d a ' a A{x, x 1 ) = B(x, x'), 



(3.7) 



where □' = ij^d^'d^; summation over a' and a in the third term is understood. By virtue of Eq. (p.6[), all time 



derivatives drop out of Eq. (3__7), which can then be written as ^(V 2 + W 2 )A(x, x 1 ) + d a a A(x,x') = B(x,x'). The 
left-hand side of Eq. ( (3.7] ) can be simplified differently. The wave operator acting on A can be taken inside the integral, 
where it acts on G flat . This returns a S- function, and the integral can be evaluated. This gives 



UA(x,x') = -2^>{x)G tlat (x,x'), n'A(x,x') = -2$(x')G Hat (x, x'), 



(3.8) 



and these expressions can be substituted into Eq. (3.7) 



B. Scalar Green's function 

We are looking for a function G(x,x') that solves |4l]] 

E[G] = -4nS 4 (x - x') 



(3.9) 



to first order in the Newtonian potential <f>. Here, Si(x — x') is the four-dimensional 5-function, and E[G] stands for 
the left-hand side of Eq. (2.7) multiplied by \J— g, with G(x,x') taking the place of ^p(x). In terms of this Green's 
function, the solution to Eq. (2.8) is 
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i/)(x)= G(x,x') t i(x')y/^g 1 d 4 x' 



(3.10) 



where g' is the metric determinant evaluated at x' . 
To find this Green's function, we write 



G(x, x') = G flat (a;, x 1 ) + G(x, x') + 0($ 2 ), (3.11) 

in which G aat (x,x') is fiat-spacetime solution given by Eq. ( |3.3| ), and G(x,x') is the first-order correction (linear in 
the Newtonian potential). [The overd ot no tation suggests that G(x,x') should be viewed as a smooth function of a 
smallness parameter e, and that Eq. (3.11) gives the first two terms of its Taylor expansion about e = 0.] By virtue 
of Eq. (|2.9|), in which we substitute □$ = Anp, we find that G(x,x') satisfies 



UG{x,x') = -^(x)d tt G tiat (x,x') + 8TT£_p(x)G ilat {x,x'). 



The solution is 



G(x,x') 



G aat (x,x")<f>(x")d t „t»G iia - t (x",x')d i x" - 2£ / G ttat (a;, x")p(x")G fiat (x" ,x') d 4 x" 



(3.12) 



(3.13) 



In the first integral, the derivative operator can be switched to dt>f because the flat-sp acet ime G reen 's function 
depends only on x" — x'. The derivatives can then be taken outside the integral. Using Eqs. ([O]) and (3^2), we obtain 
G = 2d t >t'A — 2££>. Using Eq. (3.6), we express the final result as 



G(x, x') = -2d tt 'A(x, x) - 2£B(x, x). 



(3.14) 



Equations (3.11) and (3.14) give the scalar Green's function to first order in <£>. As was claimed previously, G(x,x') 
can be expressed in terms of the generating two-point functions introduced in Sec. Ill A. For the special case £ = 0, 
our Green's function agrees with the one obtained by DeWitt and DeWitt |48| . 



C. Electromagnetic Green's function 

The electromagnetic Green's function G a p,{x,x') is a solution to jllj 

E a p,[G] = -4ir6 a ,6 4 (x - x'), 



(3.15) 



where E a p,[G] stands for the left-hand side of Eq. (2.1C) multiplied by y/—g, with G a g,(x, x 1 ) taking the place of 
A a (x). Here and below, the unprimed tensor indices refer to the point x, while the primed indices refer to x' . In 
terms of this Green's function, the solution to Eq. (2.11) is 



A a (x)= I Gyix.x')/ {x')^/^ (fx' 



We seek to determine the Green's function to first order in $. 
For this purpose, we write 



G a p, (x, x') = G flat (x, x')Sy + G a p, (x, x') + 0($ 2 ) 



(3.16) 



(3.17) 



and substitute into Eq. ( |3.15| ). Here, G tia,t (x,x') is the scalar Green's function of Eq. (3__3), and G a p, (x, x') is the 
first-order correction to the electromagnetic Green's function. We find that this must satisfy 



UG a B , = -4$<9 tt G flat (5 



flat 



(3.18) 



We have used the fact that $ does not depend on time, so that, for example, x ali ®,nf3 — V 01 ^®,^/! — such terms 
end up canceling out. We have also, on the right-hand side, dropped the distinction between (3 and j3'; this distinction 
is not necessary when dealing with a constant tensor such as x°p- 

Equation ( |3.18|) is solved by once again invoking the flat-spacetime Green's function. The solution is 



G a p/(x, x ) 



1 

An 



G^\x,x")Km{x")d A x" 



(3.19) 
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wher e RH Sfz") stands for the right-hand side of Eq. (3.1§|) evaluated at x = x" . The derivatives appearing in 



Eq. (|3.18| ) are now taken with respect to x" , When d^ii (say) is acting on G at (ic", x'), it can be replaced by — 9 ( 



and taken outside the integral. This leaves us with integrals involving G &llt (x,x"), G flat (x", x'), as well as 
and its derivati ves. Those int egrals can all be expresse d in term s of the two-point functions A{x,x') and B(x,x') 
define d in Eqs. (3.1) and (3.2), tha nks to the relations (3.4) and (3.5). The results can then be simplified by using 
Eqs. (3.6)-(3.£), as well as Eq. ( [2.3| ) for Xap- The end result is 



G> = -2d tt ,A5% + (di - d%,)A + 2t a {d t , p - d tp ,)A + 2^(0*/ - d a t ,)A + x « /3 (A$G flat - B), 



(3.20) 



where A$ = $(x) — $(a;'). Equations ( |3.17|) and ( 3.2C ) give G Q ^,(x, x'), the electromagnetic Green's function, to first 
order in 

Equation ( |3~20| ) is easier to deal with if we distinguish between its temporal and spatial components. From it we 
read off 



G\, = -A$G flat - 2d tt 'A + B, 
G\, = {d t , a -d ta ,)A, 

G\, = (d a t ,-di)A, 

G\, = 5 a b (A$G flat - 2d tt >A -B) + (d a i - d a b ,)A. 



(3.21) 



We notice that in Eqs. ( 3.20) ) and ( 3.2l| ), the part of G a p, {x,x') that involves G Rat (x,x') has support only on the 
past light cone of the point x. The remaining part, that involving the two-point functions A(x, x 1 ) and B(x, x'), have 
support inside the light cone as well; this is the "tail part" of the Green's function [^9| . Our electromagnetic Green's 
function agrees with DeWitt and DeWitt Q . 



D. Gravitational Green's function 



The trace-reversed gravitational Green's function 0°^,$, (x, x') satisfies the equation 



E a e^] = -A^ a rflMx-x% 



(3.22) 



where E a ^, s ,[Q] stands for the left-hand side of Eq. ( [2.15 ) multiplied by y/— g, but with Q al t, S / (x, x') replacing j a,3 (x). 
In terms of this Green's function, the solution to Eq. ( 2.15 ) is 

7 Q/ V) -4 / g a ?(x,x')5T"<' s '(x')^Pd 4 x'. 



(3.23) 

Our goal here is again to calculate the Green's function to first order in the Newtonian potential. The steps are virtually 
identical to what was done for the scalar and electromagnetic cases. The only difference is that the expressions are 
longer. 

At the end of the calculation we will have to "trace-reverse" the trace-reversed Green's function to obtain 
G a ^, s ,(x, x'), the Green's function directly associated with the metric perturbation 7 ai g. This is given by 



G°% (x, x 1 ) = G a % {x,x')-\ g"P (x)g, v (x)G%, s , (x, x'). 



In terms of this, the metric perturbation is recovered by integrating 



G al3 , 5 , (x, x')ST 7 ' 6 ' (x')^f^g 1 d A x'. 



To find the trace-reversed Green's function, we write 



g a ? s ,(x,x') = G^\x,x')6^5%+g a P(x,x>)+0(&) 



ict/3 



(3.24) 



(3.25) 



(3.26) 



where G aa,t (x,x') is the flat-spacetime scalar Green's function, and (x, x') the correction of order <£. By substi- 

tuting Eq. ( |3.26 ) into Eq. (3.22), we find that this must satisfy 
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- 2 



a/3 



-(flat 



- 47T 



X 7 <5 



2S^6 f3 >+4t a t^r ll s) P G 



-(flat 



(3.27) 



where we have gone through t he sa me simplification steps as in the electromagnetic case. The remaining steps are 
also similar: We integrate Eq. ( 3.27 ) with the flat- spac etime Green's function, and we simplify the resulting integrals 
as was described in the paragraph preceding Eq. ( p.20[ ) . After a rather long computation, we obtain 

g"? = -2d w A8^} + 26% \d% - d% + 2t« ( ds)t , - d s , )t ) - 2t 5) (d®, - d^)\ A + x , s (d a + d<*')(d? + d^)A 



- 2 x %(d^ + d^) (d S) + d S '))A + x aP {d, + d Y ) (d s + d s ,)A + 2 X %5%A$G 
+ (2 X %6% - V ^x,s - 25^1 - 4t"t% s )B, 



(3.28) 



where A$ = $(x) - $(x'). Equations fl3.26| ) and ( |3.28| ) give g a0 , s , (x, x'), the trace-reversed gravitational Green's 
function, to first order in $. 



From Eq. (3.28) we extract the following components: 



g" 



-2A$G 



Hat 



2d w A + B, 



t>a> ~ (dt'a- d t a')A, 

Q U a >v = (d a + d a ,)(db + d b ,)A - 35 ab B, 



g ta t , t , = (d a t ,-di)A, 



r>ta 

y t>b' 



-5\{d w A + B)+ l -(d% + 28% + di)A, 



g\, d = 5\ h {d c)t , - d c , )t )A, 



(3.29) 



g a \, t , = (d a + d a ') (d b + d b ')A - S ab B, 
6%*= 6^1 -Sty) A, 

2<5^(A$G flat ~ du>A) + S ab (d c + d c ,) (d d + d d ,)A - 25 {a [c (d b] d) + 2d b) dl) + 8% )A 
+ 5 cd (d a + d a ') (d b + d b ')A - 6 ab 6 cd B. 



r>ab 

y c'd> 



To "trace-reverse" the Green's function according to Eq. (3.24) is a straightforward operation, but it is important to 
keep in mind that the metric involved in this computation is not rf 1 ^ , but rj a ^ + 2$x a/3 . We find that the gravitational 
Green's function is given by 



G af3 Y5 ,(x,x>) = (6%S f3 l - \r ] a %, 5 )^\x,x') +G^, 5 ,{x,x') + 0^)\ 



where 



The complete listing of components is 



G u t , t , = -(A$G nat + dtfA), 
G u t , a , = (d t , a -d ta ,)A, 

G u a , v = -5 ab [{$)G aat + d tt ,A + 2B] + {d a + d a ,)(d b + dv)A, 



(3.30) 



G a % s ,{x,x>) = g af) YS ,(x,x>) - -v^G^rpM - 2{t a t%, s , - r%>U')Hx)G^{x,x>). (3.31) 



G ta vt , = (d a t ,-di)A, 



Gta 
t'b' 



-6%(d tt 'A + B) + -{d\ + 2d\ + d\,)A, 



G\, c , = 8 a {h (d c)t ,-d c , )t )A, 



(3.32) 
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G a \, t , = <5 afc [($)G flat - d w A] + (d a + d a ') (d b + dv)A, 

G ab dd , = (25^6% - 5 ab 6 cd ) (A$G flat - fa A) + S ab (d c + d c ,)(d d + d#)A 

- 26 {a (c {d b) d) + 2d h) dl) + d b 'J, } )A + 5 cd (d a + d a ') (d b + d b ')A - 25 ab S cd B. 



Here we use the notation ($) = + as well as A$ = — $(x'). It should be noted that in Eq. (3.32), the 

part of G a P, g ,(x, x') that involves G Rllt (x, x') has support only on the past light cone of the point x. The remaining, 
tail part has support inside the light cone as well. We believe that these results for the gravitational Green's function 
are new. 



IV. EVALUATION OF TWO-POINT FUNCTIONS 



The scalar, electromagnetic, and gravitational Green's functions can all be expressed in terms of A(x,x') and 
B(x,x') 1 the two-point functions introduced in Sec. Ill A. Our task in this section is to evaluate these two-point 
functions. We will keep $>(x) completely generic, demanding only that far from the matter distribution, it becomes 
equal to —M/r, where r — \x\ and M is the total mass. This distinguishes our work from the earlier work of DeWitt 
and DeWitt p8| , who considered only the special case $ = —M/r. A consequence of keeping $ generic is that 
we will obtain only partial information regarding the two-point functions. This, however, will be sufficient for the 
computation of the self-forces presented in Sec. V. The results specifically required for this computation are derived 
in this section. Additional results are presented in Appendix A, which gives an extended discussion of the two-point 
function A(x,x'). Complete expressions for the two-point functions are derived in Appendix B for the special case 
<t> = —M/r; these reproduce the earlier results of DeWitt and DeWitt. 



A. Ellipsoidal coordinates 



When Eq. (3J3) is substituted into Eqs. (3.1) and (3.2), we find that the dt" integration can be carried out imme- 
diately, and we obtain 



A(x,x') = 



1 

27 



$(cc") 



IX — X 



— 6 (At - \x - x"\ - \x" - x'Q d 



V 



(4.1) 



and 



B(x,x') = J ■ 



P (x") 



— 6(At - \x - x"\ - \x" - x'\) d 



V, 



(4.2) 



where At = t — t' . These are the integrals we will attempt to evaluate. The method of calculation presented below 
comes from DeWitt and D eW itt E^L but this particular implementation was suggested to us by Alan Wiseman p5| . 

The (5-function in Eqs. (4.1) and (4.2) enforces the condition \x — x"\ + \x" — x'\ = At, which defines the closed, 
two-dimensional surface S formed by the intersection of x's past light cone with x n s future light cone. This surface is 
an ell ipso id of r evo lution centered at Xq = | (x + x'), of semi-major axis \At and ellipticity ||x — x'\. To integrate 
Eqs. (4T) and (O), we will adopt ellipsoidal coordinates adapted to this geometry. 

Wc first summarize our notation: 



At = t- t', 



Xq = ~(X + X), 



R = x — x' , R = \x — x'\ 



n=*, e=\R (4.3) 



The diagram of Fig. 1 illustrates the situation. 

For convenience we orient our coordinate axes such that the unit vector ft points in the same direction as the z, or 
third, axis; there is no loss of generality involved in this choice. We then represent the vector x" in terms of ellipsoidal 
coordinates is, 9, 6). The coordinate transformation is 



x" = Xq + T](s, 



(4.4) 
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FIG. 1. The diagram shows the matter distribution located near the origin of the Cartesian coordinate system, the points 
x, x' , and xq, to which the ellipsoidal coordinates (s,8,4>) are attached. The diagram also shows one representative ellipsoid 
of constant s, and a point x" — xq + rj(s, 8, <f>) lying on this ellipsoid. 



where 



n 



sin V cos < 



sin V sin < 



V 



(4.5) 



The parameter e 



ji? is the ellipticity of the new coordinate system. It is easy to check that in these coordinates, 



s — e cos 9 and Ice' 



x \ = s 



e cos 9. The Jacobian of the transformation ( [4.4|) is (s 



\x — x i 

With these results, we find that Eq. JO) reduces to = (2tt)^ 1 J ^>(x Q + rj)8(At 

sin 9 d9d<j) is the element of solid angle. The integration over s is immediate, and we obtain 

A(x,x') = ±- { ^(x + V (9^)) dSl, 



I sin f 



2s) dsdfl, where dfl 



(4.6) 



is now given by Eq. (4.5) with s set equal to iAt. Thus, A(x,x') is the average of $ over 



where the vector r)(9 

the ellipsoid s = -kAt, which is just the two-surface S introduced previously. Similarly, 



B(x,x') 



p(x + r)(6,4>))dn 



(4.7) 



is the average of 2irp over the same ellipsoid. It is important to remember that in Eqs. (4.5)— (4.7), the third coordinate 
axis has been aligned with the vector R = x — x'. 

It should be noted that A(x, x') and B(x, x') are nonzero if and only if x is in the causal future of x'. To properly 
reflect this, the integrals of Eqs. (4.6) and (4.7) should be multiplied by the step function 9(At — R). To keep the 
notation simple, however, we choose to leave this factor implicit in our expressions; for the rest of this section, At 
will be understood to always be larger than R. 



B. A(x,x') near coincidence 



The integral of Eq. (4.6) can readily be evaluated, for an arbitrary Newtonian potential if the points x and x' 
are close together in spacetime. This is the near-coincidence approximation. Because these points must be timclike 
related for A(x,x') to be nonvanishing, we are looking at evaluating Eq. ( p~6| ) in a limit where R and At are both 
small, but the ratio At/R is maintained larger than unity. In this limit, the vector r/(9, <fi) defined in Eq. (4.5) is small, 
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and $ can be expanded in Taylor series about the point Xq; this is the basis for our approximation. For concreteness, 
we will assume that Xq lies outside the ma tter distribution, so that p{xq) = 0. 
Expanding <f> in powers of r\ in Eq. (4.6) gives 



4ttA = <P(x ) [dn + $, a (xo) j V dfl+^, ab (x ) J ?yV dQ 



(4.8) 



The first integral, with a factor of (47r) 1 in front, gives unity. The average of r) a on the sphere is zero, and the average 
of rj a r] b is the tensor q ab = |(s 2 — e 2 )S ab + ^e 2 h a h b , where s = |At and e = The part of q ab that is proportional 
to S ab multiplies &, a b(xo) and returns something proportional to V 2 <£> evaluated at Xq; because this point is outside 
the matter, this is zero. We are left with 



1 



A = S(s ) + 24 R 2 ^,ab(xo)n a h b + 



(4.9) 



Notice that the right-hand side of this equation no longer depends on At: the two-point function A(x,x') is time- 
independent near coincidence. 

It is convenient to shift the reference point from Xq to x. To order R 2 , this means re-expressing $>(xq) as 



<f>(x - \Rh) = - - m, a {x)h a + - R 2 $ M (x)h a n b + ■ 



Substituting this into Eq. (4.9), we arrive at 



A(x,x') = $(x) - i $Ax)R a + \ <$>, ab {x)R a R b + 0(R 3 ). 
2 6 



(4.10) 



(4.11) 



This is our final expression for the two-point function in the near-coincidence approximation. We recall that R = x—x'. 

In the sequel we will need expres sions for various derivatives of A(x,x') evaluated at coincidence, x — x' . It is 
straightforward to differentiate Eq. ( 4.11 ) with respect to either x or x', to obtain 



lim d t A = lim d v A = 



(4.12) 



and 



lim d ab A = 2 lim d ab >A = lim d a >vA = - $ iab (x). 



(4.13) 



These results hold for arbitrary Newtonian potentials, provided that the point Xq lies outside the matter distribution. 



C. A(x,x') for long delays 

We can also evaluate A(x,x') for arbitrary potentials if we assume that x and x' are widely separated in space- 
time. More precisely, we now take At to be extremely large, and in particular, much larger than R = \x — x'\, 
r = \x\, and r' — \x'\. This is the long-delay approximation. In this limit, s is much larger than e, and 



77 ~ (s sinO cos (f>, s sin 9 sin <p, s cos 9). Because s is very large, we will neglect Xq in front of 77 in Eq. (4.6), so 
that A ~ (Att)^ 1 J &(r]) dCl. Since 77 is located well outside the matter distribution, we can approximate the Newto- 
nian potential by —M/\ij\, its first term in a multipole expansion. (Here, M = J pd 3 x is the total mass of the matter 
distribution.) Because \r/\ ~ s = \At in the long-delay limit, we arrive at 

A(x,a/)^~. (4.14) 



We recall that At = t — t'. Equation (4.14) shows that A(x, x 1 ) and all its derivatives vanish in the limit t' — > —00. 
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D. B(x,x') at large distances 



We will evaluate B(x, x') in the case w here x and x' are both located well outside the matter distribution. For this 
calculation it is easier to go back to Eq. (4.2), in which we assume that r = \x\ and r' = \x'\ are both very large. (A 
more precise criterion will be given below.) These conditions will allow us to develop a multipole expansion for this 
two-point fun ctio n. 

Equations (4^2) and (4/7) indicate that B(x,x') is the average of the mass density p over the ellipsoid s = \At. 
The diagram of Fig. 1 makes it clear that unless At falls within the appropriate interval, the ellipsoid will fail to 
intersect the matter distribution. Thus, B(x, x 1 ) is nonzero only when At lies within this interval, which we might call 
the intersection window. When r and r' are both large, the ellipsoid is also very large, and the intersection window 
becomes comparatively short. In this limit, B(x,x') is "on" for a very short time, and its behavior suggests that of a 
5-function. This expectation is b orne out by an explicit calculation. 

In the denominator of Eq. (4.2), we approximate \x — x"\ by r and \x" — x'\ by r', taking \x"\ to be much smaller 
than both r and r' . In the <5-function we go to the next order of approximation and write \x — x"\ ~ r — x ■ x"/r 
and \x" — x'\ ~ r' — x' ■ x"/r'. If we introduce the vector 



x x 
m = — | — - 

r r 



(4.15) 



then the (5-function becomes S(u 



m 

5(u) + 5'{u){m ■ x") + \8"(u)(m ■ x"f + 
Substituting all this inside the integral of Eq. (4.2), we obtain 



'), where u = At — r — r' . This we expand in powers of x" and express as 
wh ere a prime denotes differentiation with respect to the argument. 



B 



1 

rr' 



S(u) I p(x")d 3 x" + S'(u)m a I p{x")x a " d 3 x" + ^-8"{u)m a m b I p(x")x a " x b " d 3 x" 



(4.16) 



The first integral gives the total mass M. The second integral gives the dipole moment of the mass distribution; this 
vanishes if, as we assume, the origin of the coordinate system is attached to the center of mass. The third integral is 
Q ab , the quadrupole moment of the mass distribution. We have therefore obtained 



1 



B(x,x') ~ — 



1 



M5{u) + -Q ab m a m b 5"(u) + 



(4.17) 



where, we recall, u = At- 



This equation displays the first two terms in the multipole expansion of the two-point 



function B(x,x'). Notice that this expansion is analogous to a multipole expansion of a radiative field in the wave 
zone; this is not a near-zone expansion of a quasi-stationary field. 
In the sequel we will approximate B{x,x') by 



M 



B(x,x') = — 5 (At 



0. 



(4.18) 



its leading term in the multipole expansion. In this approximation, B(x,x') is nonzero only when At = r + r' . The 
time delay corresponds to a signal propagating at the speed of light from x' to the center of the mass distribution, 
and then on to s. The approximation has made the i nterse ction wi ndow infinitely short. 

To estimate the error involved in going fro m Eq. ( 4.17 ) to Eq. ( 4.18 ), we integrate B(x,x') over a smooth test 
function f(t). Setting t' = for simplicity, Eq. (4.17) gives 



Bf(t)dt~^jf(r- 



1 



M f 



(4.19) 



whe re Q = ^Q a t > m a m b . The factor in front the large brackets is what would be returned by the approximation of 
Eq. J4.16| ). We wish to show that in typical situations, the second term inside the brackets is much smaller than unity. 
Let I be the characteristic size of the matter distribution, and let r be the characteristic time over which the function 
f(t) changes. Then Q/M is of order £ 2 , and /"// is of order 1/r 2 . The correction term is therefore of order (£/t) 2 . 
In a typical situation, r would be associated with the orbital time scale of an object located at x or x'. Assuming 
for concreteness that r is the shortest distance, we have r 2 ~ r 3 /M. With this, we find that the correc tion term in 
Eq. ( 4.19 ) is of the order of (£/r) 2 (M/r). So long as r and r' are both much larger than I and M, Eq. ( 4.1S ) makes 



an excellent approximation. 
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V. SELF-FORCES 



A. Scalar self-force 



In the absence of an external force, a point particle with scalar charge q experiences a self-force [|46| 

r(r) = q 2 [g af3 + u a (r)u^r)] T G i( g(r, r') dr' , (5.1) 



so that its equations of motion are mv^.puP = f a . The particle's world line is represented by the relations z a (r), 
where r is proper time, and u a (r) — dz a /dr. The self-force is evaluated at the current position of the particle, and 
it involves an integration over its past history. Inside the integral, the retarded Green's function is first written as 
G(x, z(r')), in terms of an arbitrary field point x and the past position z(t') of the particle. The Green's function is 
then differentiated with respect to x@ , and th e re sult is evaluated at x — z(r), the particle's current position. This 



is what we have denoted by G i/ g(r, r') in Eq. (5J). The upper limit of integration is t~ = r — e, where e is a small, 
positive number that is taken to zero at the end of the calculation. The purpose of this cutoff is to remove the singular, 
(5-function part of the Green's function which has support on the past light cone of the point z(t). What survives is 
the "tail part" of the Green's function, which has support inside the light cone. 

This expression for the self-force was derived by Quinn |46|| , who assumed that the scalar field was minimally 
coupled to curvature — he considered the special case £ = 0. There is no difficulty, however, in extending Quinn's 
result to arbitrary couplings to curvature. So long as the particle moves in a region of spacetime in which the Ricci 
scalar is zero, his results carry over directly. The only effect that £ has on the motion is through the Green's function, 
which is sensitive to a nonzero Ricci scalar somewhere in spacetime. 

We will calculate f a to first order in the Newtonian potential $ in a weak-field approximation; the particle therefore 
moves in the weakly curved spacetime of Sec. II. We assume that the particle is gravitationally bound to the matter 
distribution, and infer from the virial theorem that v 2 , the square of the particle's velocity, is of the same order of 
magnitude as $. We therefore add a slow-motion assumption to our weak-field assumption. Accordingly, we will 
neglect in all expressions terms involving <i> 2 , $u 2 , and u 4 . Within this approximation, we have that the spatial 
components of the self-force are given by 

f a =q 2 I (G ia + v a G it )dt', (5.2) 



where v = dz(t)/dt is the current velocity of the particle, which does not depend on t'; the past velocity will be 
denoted v' = dz(t')/dt', and this does depend on t' . We have made t' the new variable of integration, and neglected 
all 0{v 2 , vv' , v' 2 ) terms inside the integral, which through the Green's function is already linear in the Newtonian 
potential. This fact allows us to freely alter the position of all latin indices, implicitly using S a b, the metric of fiat 
space. And since we will substitute only the tail part of the Green's function inside the integral, we have made t the 
upper limit of integration (instead of the redundant t — e). 



From Eq. (3.14) we find that the tail part of the scalar Green's function is given by 

G{x, z') = -2d tt 'A(x, z') - 2£B(x, z'). (5.3) 
Here z' = (t' , z(t')), and A(x, z'), B(x, z') are the two-point functions introduc ed i n Sec. Ill A and computed in Sec. IV; 



we will set x — z = (i, z(t)) after taking the derivatives. Substituting Eq. (5^) into Eq. ( |5.2| ) gives f a = f% + f%, 
where 



r t 

-2 



Ja = - V / (A,tt> a + v a A M , t ) dt' (5.4) 

and 

f a B - -2£<? 2 / ill, - r'-B.ndl I -V.-.l 
We will first work on f\, and then move on to f%. 



The first term inside the integral of Eq. (5.4) can be expressed as d t >A ta ~ dA ta /dt' — A^yv , allowing part of the 



integrati on to be carried out. Similarly, we write the second term inside the integral as v a dt'A^t = v a dA tU /dt' +0(vv'), 



and Eq. (5.4) becomes 
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n = 



A, ta + v a A t , 



2q 2 



A M b>v b ' dt'. 



(5.6) 



We now use Eq. (3.6) and replace A tta b' by —d t 'A. a y — —dA^/dt' + 0(v'). After integrating by parts, we obtain 



A, ta +v a A. tt + A. ab ,v b ' 



A n wa h ' dt' 



(5.7) 



where a' = dv' /dt' is the past acceleration of the particle. Because the particle is not subjected to any force except 
for g ravi ty, its acceleration is linear in the Newtonian potential Since A(x, z') is already linear in the integral of 
Eq. ( |5.7| ) is 0(4> 2 ) and it can be neglec ted. The first part of the self- force can therefore be computed without having 
to evaluate a single integral! And Eq. (5.7) shows that the self-force depends on the conditions at the current time, 
and in the infinite past. At this stage, the nonlocal character of the force has (almost) disappeared. 

In Sec. IV C we learned that as t' approaches — oo, A(x,z') behaves as 2M/t'. This shows that in Eq. (5.7), there 



are actually no contributions from the infinite past. Any nonvanishing contribution must therefore come from the 
current time, t' = t, at which the points z(t') and z(t) are coincident. Equation (4.12) tells us that the time derivative 
of A(x, z') is zero at coincidence. Setting v b = v b , we find that the sole surviving contribution to f% is —2q 2 A, a b> v b , i n 



which the derivatives of the two-point function must be evaluated at coincidence. This was worked out in Eq. (4.13) 
and we arrive at 



n 



We may express this in the alternative form 



fA 



3 9 dt' 



(5.8) 



(5.9) 



where g = — V$ is the Newtonian gravitational field. In components, the right-hand side of Eq. ( |5.9| ) involves 
dg a /dt — —d<fr :a /dt = —$ !a bV b , where we have used the fact that the Newtonian potential does not explicitly depend 
on time. 

We now turn to the second part of the self-force, given by Eq. ( |5.5| ) in terms of the two-point function B(x,z'). 
According to Eq. ( 4.18] ), this is 



M 



B(x,z') = —S(u) 



(5.10) 



where u — t — t' — r — r' , r = \x\, and r' = |z(t')|- As before, the field point x = (t,x) will be identified with 



z = (t, z(t)) after differentiation. It is important to keep in mind that in Eq. ( 5.1C ), r' depends on t' . We will use the 
notation v' r — dr' /dt'. 

We first show that in Eq. ( [5.5| ), the integral involving the time derivative of B(x 7 z') is negligible. According to 
Eq. ( |5~T(i| ), B tt = M6'(u)/(rr'), where a prime indicates differentiation with respect to u. Changing the variable of 
integration from t' to u, using du = —(1 + v' r ) dt' , we have that / B^dt' = (M/r) /(r') _1 [l + 0(v')]S'(u) du. Here 
the lower limit of integration is — (r + r'), the upper limit is oo, and we neglect terms of order v' because the integral 
comes with a factor of v a in front. Integration yields Mv' r /(rr 12 ) evaluated at u — 0. We conclude that the second 
term of Eq. (5.5) is of or der ( q/r) 2 (M/r)v 2 and therefore negligible. 

Differentiation of Eq. (5.1C) with respect to the spatial variables x a yields 



B, a = - 



Mr a 



6(u) + rd'(u) 



(5.11) 



where a prime still indicates differentiation with respect to u. Changing variables of integration from t' to u, we find 
that the first integral of Eq. fl5.5|) becomes 



B n dt' 



Mr, a 
Mr, a 
Mr, a 



— (r+r 1 ) 
1 



5{u) + rd'(u) 
r'(l+v' r ) 
d 

— r- 



du 
1 



r'(l + v' r ) du r'(l + v' r ) 
1- (l + r/r')v' r + 0(v' 2 ,a r 



u=0 



u=0 



(5.12) 
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where, as indicated, we neglect terms that are of the second order in the velocity, and terms that involve the particle's 
acceleration. 

To evaluate the right-hand side of Eq. ( |5.12 ) is made difficult by the fact that the relation u = t — t' — r — r'(t') 
cannot be inverted for t' , because the function r'(t') is not known. For the purposes of this calculation, however, it 
is sufficient to give an approximate inversion. We let the radial part of the particle's trajectory be described by the 

relation r = R(t), so that r' = R(t') = R(t - At). This we approximate by R(t) - R(t)At + \R(t)At 2 H , where 

overdots indicate differentiation with respect to t. Neglecting the acceleration term, we have r' = r — v r At + 0(a), 
where v r = R(t). Similarly, we obtain v' r = R(t') = v r + 0(a). When u = 0, we fin d tha t At = 2r j (1 + v r ) + 0(a) = 



2r(l — v r ) + 0(v 2 ,a), giving r' = r(l — 2v r ) + 0(v 2 ,a). Substituting this into Eq. ( 5.12) ), we find 



B a dt' = 



Mr„ 



l + 0(v 2 ,a) 



(5.13) 



Finally, going back to Eq. (p.5|), we arrive at 



f B = 2^r, 



(5.14) 



where f — x/r is a unit vector pointing in the radial direction; its components are equal to x a /r 



The total self-force acting on a point particle with scalar charge q is given by the sum of Eqs. (5.9) and (5.14). This 



is 



,dg 

dt 



(5.15) 



where, we recall, g = — is the Newtonian gravitational field. The equations of moti on for the particle are then 
d 2 z/dt 2 = g+ f /m, where m is its mass. The first term on the right-hand side of Eq. ( |5.15| ) is a conservative correction 
to the local gravity. This is a repulsive force if £ > 0, an attractive force if f < 0, and this part of the self- force 
disappears altogether if the scalar field is minimally coupled to curvature. The second term is nonconservative; this is 
the well-known expression for the radiation-reaction force experienced by a scalar charge moving under the influence 
of an external force mg. 



B. Electromagnetic self- force 

The steps required to calculate the electromagnetic self-force are very similar to those presented in the preceding 
subsection. Our discussion here will therefore be brief. We employ the same notation, and work within the same set 
of approximations. 

Our starting point is the DeWitt-Brehme expression for the self-force experienced by a point particle with an electric 
charge e El|: 



r = -e 2 



{G^.p-G^uW dr', 



(5.16) 



where G a g,(T, r') is the retarded Green's function of the electromagnetic field, u a = dz a (r)/dr is the current four- 
velocity of the charged particle, and u a = dz a (T')/ 'dr' its past four- velocity. With the approximations introduced in 
Sec. V A, we find that the spatial components of the self-force are given by 



r 



(G a t',t — Gtt',a) + (Gat',b — Gbt',a)v b + (Gafc'.t ~ G t b\a)v b 



dt'' 



(5.17) 



only the tail part of the Green's function appears under the integral. Using Eq. (3.21), we find that the self-force can 
again be written as a sum of two parts, / = /a + /s, where 



n = e 2 



(Ajt'a + Atta') + (Aja'b - A Mb ,)v b + (2S ab A, tt 't ~ Ata'b + 2A Mb , - Aj, ab )v b 



dt' 



(5.18) 



and 
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f% = 



j (B, a - B t v a ) dt'. 



(5.19) 



The evaluation of fs proceeds as in Sec. V A, using Eq. ( 5.13 ) and the fact that the second term does not contribute 
to the force at this level of approximation. The result is /§ = e 2 (M / r 3 )r M . The evaluation of /a also involves familiar 
steps, which we now describe. 



The first group of terms inside the integral of Eq. (5.1 8|) are expressed as 



d 



d t >(A M - A tta >) = -jpiAta ~ Ata>) ~ {A, ta b> ~ A ta , b ,) V h 



(5.20) 



The first term is immediately integrated, and by virtue of Eq. ( |4.12| ), this gives no contribution to the self- force. The 
remaining terms are absorbed into the third group of terms in Eq. ( 5.18 ), those which are proportional to v b . In the 
second group of terms, the quantity in front of v b is expressed as 



d v (A^ b - A ta b>) = --Tjj(A,a>b - A, ab ,) + 0(v'). 



(5.21) 



This is also immediately integrated. Equation ( |4.13| ) shows that this also gives no contribution to the self-force. We 
are therefore left with the group of terms proportional to v b . This takes the form of (dC a b/dt')v b , where 



Cab — 25 a bA : tt — A^b ~ A^ a y + A^ a 'b — Ay 



(5.22) 



The time derivative is then transferred to v b by integration by parts. This produces an integral that can be neglected, 
and boundary terms at t' = t (the boundary terms at t' — — oo all vanish). Those give C a bV b evaluated at coincidence. 
Using Eqs. (f4.12|) and (f4.13|) we find that Eq. Q5.18|) becomes f% 



§e 2 $ ah? A 



We have found that the self-force acting on a point particle with electric charge e is given by 



2 M . 2 2 „ ; , 



dg 

dt ' 



(5.23) 



where r = x/r and g = — V<I> is the Newtonian gravitational field. The first term on the right-hand side of Eq. ( 5.23| ) 
represents a repulsive, conservative force. The second term gives the well-known expression for the radiation-reaction 
force experienced by an electric charge moving under the influence of an external force mg |Jl3"| , fl4f . 

Equation (5.23) implies that the force required to keep an electric charge static in the gravitational field of an isolated 
object of mass M must be smaller than mg by e 2 M/r 3 . Our derivation indicates that this correction is caused by 
an interaction between the electromagnetic field and the matter distribution; this interpretation is suggested by the 
fact that this contribution to the self-force can be traced back to the two-point function B(x,x'), which is directly 
associated with p, the mass density. While this is a valid interpretation, it is interesting to note that the correction 
persists even in the total absence of matter, although this cannot be revealed in a weak-field approach. The force 
required to hold an electric charge in place in the exact field of a Schwarzschild black hole was calculated by Smith 
and Will ]49[ |. In the weak-field limit, these authors recover the conservative term in Eq. ( 5.23| ). In this case, clearly, 
the effect has nothing to do with the presence of matter. It must instead be attributed to an interaction between 
the electromagnetic field and the hole's event horizon. That the conservative self-force should be so insensitive to the 
source of the gravitational field is quite interesting. This robustness, however, doesn't seem to apply to scalar charges: 
The extra force required to keep a scalar charge in place in a Schwarzschild spacetime is always zero, irrespective of 
the coupling to curvature |^,^5|-|8|. This is quite different from what was obtained in Sec. V A, where we found that 
in the presence of matter, there is a nonzero self-force for nonminimal coupling. 



C. Gravitational self- force 



The self-force experienced by a point particle of mass m moving in a curved spacetime is given by p6| , p7| 

/ Q - -2m 2 f (2G V,, ;7 " G Millv ; a + u a Gp ltllv ,,su s ) u^u^u^'u v ' dr' , (5.24) 



where G a ^ , v , (t, t') is the gravitational Green's function, u a the current four-velocity of the particle, and u 11 its past 
four- velocity. This result was derived for the first time by Mino, Sasaki, and Tanaka [Eq], although an incomplete 
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attempt was made earlier by Morette-DeWitt and Ging |5^] . The gravitational self-force was later recovered by Quinn 
and Wald p7fl , who gave a more transparent derivation. Their expression is technically wrong, however, because it 
incorrectly puts G „i v i, the trace-reversed Green's function (see Sec. Ill D), in place of G a ^ , v , in Eq. (5.24); this 
slight oversight hardly diminishes the value of their paper. 

Our expression for the gravitational self-force differ s from Mino et al. p6[ in two ways. First, we use a different 



normalization for the Green's function, and Eq. (5.24) is valid in the normalization imposed by Eq. (3.25). Second, 
and more importantly, both Mino et al. and Quinn & Wald |B7| assume that the background metric g a p is a 
solution to the Einstein field equations in vacuum, a condition that is incompatible with the weak-field appr o ximat ion 
adopted in this paper. [Please refer back to the discussion presented in the paragraphs before Eq. ( 1 . 13| ) .] As we 
have already indicated, the extension of the Mino-Sasaki-Tanaka-Quinn-Wald equations of motion to spacetimes that 
contain matter is not entirely trivial. In this subsection we consider only the s traightforward modifications to the 
Green's function that come fro m th e presence of Ricci-tensor terms in Eq. (2.13). In other words, we still define the 
gravitational self- force by Eq. ( 5.24 ), but we substitute in it the Green's function computed in Sec. Ill D to account 
for the presence of matter. In Sec. VI we shall prove that this substitution is appropriate, and derive additional 
(matter-mediated) corrections to the equations of motion. 

Under the weak-field, slow-motion approximation, we find that the spatial components of the self-force are given 

by 



f a = -2m 2 



{2G a tt't',t — Gttt'V ,o) + (8abGttt't',t + %GaH't',t + 2G a tt't',b ~ 2Gbtt't' ,a) vb 



+ 2(2G a tb't',t — Gttb't',a)v b dt' 



With the tail part of the Green's function extracted from Eq. (3.32), we have f = fA + /s, where 



(5.25) 



(2A <tta > - A <tVa ) + (-3S ab A. tt , t + 2A Mb + 4A <ta , b + 2A M , b ,)v b 



+ 2{28 ab A M 't - 2^, tQ 'b - A M >v + A Mb >)v b 



dt' 



(5.26) 



and 



B t v a dt'. 



(5.27) 



Referring back to the discussion preceding Eq. ( 5.11 ), we see that fs is of the order of (mM/r 2 )(m/r)v 2 , which is a 
correction of order relative to Newtonian gravity. This gives a second post-Newtonian (2pn) correction [g2| to 
the Newtonian equations of motion, and according to our approximation rules, this must be neglected. 

The evaluation of f& is lengthier, but the steps are familiar. The first group of terms inside the integral of Eq. ( 5.26| ) 
are 



d v (2A M , +A. ta ) = - — {2A M , + A tta ) + (2A )to /y + A ttab ,)v h 



(5.28) 



The total derivatives contributes nothing to the self-force, and the remaining terms are absorbed into the third group. 
In the second group, the quantity in front of v b is 



dt 



7 (-3<UA«' + 2A ,ab + ^A :a , b + 2A >a , v ) + 0(v'), 



(5.29) 



and this group contributes 4m 2 $ iCl bU b to the self-force. Finally, the third group of terms is 

2v b ' -^(2SabA, tt + 2A, a , b - ^A, ab ,) + 0(v' 2 ), 



(5.30) 



and after integrating by parts, we find that this contributes — im 2 $ ia (,u b to the self-force. Summing the contributions 
gives f% = ^-m 2 $ afc iA 



The gravitational self-force acting on a point particle of mass m is therefore 



(5.31) 
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where g = — V<I>. This gives a correction of 1.5pn order to the equations of motion, and appears to represent radiation 
reaction. According to the naive equations of motion d 2 z/dt 2 = g + f /m , this term does work on the particle at an 
average rate dW/dt = +^-m 2 \g\ 2 . We see that the self-force of Eq. (5.31) seem s to bring energy to the particle. 

There are two major problems of interpretation associated with Eq. ( 5.31 ). First, the self- force seems to give 
rise to radiation antidumping instead of radiation damping Second, radiation reaction seems to occur at 1.5PN 
order, while it is known that the effect should appear only at 2.5pn order in a post-Newtonian expansion of the 
relativistic equations of motion J3|-p|,pl|. In the next section we shall see that these problems disappear once we 
properly incorporate additional (matter-mediated) corrections to the equations of motion. 



VI. MASSIVE PARTICLE IN A SPACETIME CONTAINING MATTER 



Our task in this section is to produce a proper extension of the Mino-Sasaki-Tanaka-Quinn-Wald equations of 
motion [|26|,|27|] to spacetimes that contain matter. Fi rst (Sees. VI A and B), we establish that in this context, the 



gravitat ional self- force / s " lf is still given by Eq. (L3), but that the retarded Green's function must be a solution 
to Eq. (3.22), which includes Ricci-tensor terms generated by the matter. That this does indeed give the correct 
expression for the self-force was presented as an assumption in the preceding sections of this paper. Second (Sec. VI 
C), we derive the existence of an additional term in the equations of motion, which now take the form 



;/3< 



,P - 



PCX I POL 

/self "T Jmr 



(6.1) 



where /^ m is the matter-mediated force alluded to previously. This force originates from the change in the background 
metric associated with the motion of the mass M around the system's center of mass; this motion is caused by the 
gravitational action of the particle on the central mass. Third (Sees. VI D to H), we compute the matter-mediated 
force, and show that it contains a radiation-damping term that precisely cancels the antidamping term in the self-force. 
Such a cancellation was noticed a long time ago by Carmeli |5q ], but in the context of a very different formulation 
of the equ atio ns of motion. We will see (Sec. VI I) that within the weak-field, slow-motion approximation considered 
here, Eq. (6.1) reduces to the appropriate limit of the standard post-Newtonian equations of motion &-||,0|. 



A. Background spacetime and test particle 

Let us consider the background spacetime first. We are given a metric g a p that satisfies the Einstein field equations 
in the presence of matter. We assume that the matter distribution is bounded and describes an isolated object of mass 
M; throughout this section we will refer to this object as "the star" . We write the field equations as G a ^[g] = 8nT a P[g], 
and indicate explicitly that the stress-energy tensor is a functional of the metric. For the moment we take g a $ to be 
an exact solution to the field equations; the weak-field approximation will be incorporated at a later stage. And for 
the moment we avoid adopting a specific phenomenology for the stress-energy tensor; the star can always be thought 
of as a ball of perfect fluid with a specified equation of state. 

In this spacetime we insert a point particle of mass m. We assume that this particle moves on a world line z a (r) in 
a region of spacetime that is empty of matter; the parameter r is the particle's proper time. We write the particle's 
stress-energy tensor as 

r al3 {x) =m J u a vP {~g)- 1/2 5 4 {x - z{t)) dr, (6.2) 

where u a (r) — dz a /dr is the four- velocity and the <5-function is a four-dimensional distribution with support on 
the world line. It is easy to show that conservation of the stress-energy tensor in the background spacetime implies 
geodesic motion: 

r a0 p = mj u a ^u li {-g)- x l 2 8±{x - z(t)) dr, (6.3) 

so that t q/3 . /3 = leads to u Q /3 u /3 = 0, the geodesic equation. Here, a semicolon indicates covariant differentiation in 
the background metric. We will keep this notation below: a semicolon will always refer to the metric g a fi, which will 
always be used to raise and lower Greek indices. 

We see that the particle's stress-energy tensor involves the background metric, through the definition of proper time 
and the factor (— g)^ 1 ^ 2 ■ If we think of expressing the metric as a formal expansion in powers of M, we can similarly 
expand the particle's stress-energy tensor as 
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T afi = 0{m) + 0(mM) + 0{mM 2 ) + ■■■. (6.4) 

To keep track of the powers of m and M appearing in various quantities will ex trem ely important in the following 
discussion. For the weak-field application considered here, we shall truncate Eq. ( |3.4| ) after the first two terms. 

B. Self- force 

The picture described thus far is that of test mass moving on a geodesic of the background spacetime: the particle 
has not yet been allowed to alter the spacetime's geometry. We now incorporate this effect by inserting r a " on the 
right-hand side of the Einstein field equations. Working linearly in the small mass m, we find that this modifies the 
metric by a term 7 Q /3, so that the total metric is now g a p + 7 Q /3. If we impose the Lorentz-gauge condition on 7 Q ^, 



as we have done in Sec. II B, we find that the metric perturbation is given by Eq. (3.25) 



7°"V) = 4 / G a %(x,x')^^'(x')^d 4 x', (6.5) 



where G a/3 , v , (x, x') is the gravitational Green's function whose "trace-reversed" counterpart is a solution to Eq. (3.22) 



— Green's equation for a spacetime containing matter. If we again think of expressing the metric as a formal expansion 
in powers of M, we can schematically write 

G a ^, u , = 0(1) + O(M) + 0(M 2 ) + ■■■, (6.6) 



which corresponds to Eq. ( 3.30| ). Expanding Eq. ( |6.5| ) in a similar way, using Eqs. (6.4) and ( |6.6| ), we obtain 



7ct/3 = 0(m) + O(mM) + 0(mM z ) + ■■■. (6.7) 

This expansion omits terms of order m 2 (and higher) that would come from the nonlinearities of the Einstein field 
equations. We shall systematically omit such terms in future expansions, and consistently work to first order in the 
small mass m. 

The gravitational self-force corresponds to the metric perturbation j a p acting on the particle. If this were the only 
metric perturbation to be considered — and we will show below that it is not — this action could be described by 
the statement that the particle now moves on a geodesic of the full metric g a p + J a f3- Working again to linear order 
in the perturbation, we find that the new equations of motion read 

<X = ~k-f a0 + 7V " + ^7«aXK^ (6-8) 



and the right-hand side can be identified with f^/m, the gravitational self-acceleration. Based on Eq. (6/7), we 
would naivel y co nclude that f^/rn = 0(m) + 0(mM) + 0(mM 2 ) + ■ ■■. 

Equation ( |6.8| ), however, is not valid as it stands. First, we have already indicated that we must insert a matter- 
mediated force on the righ t-hand side. Second, we know that j a p is singular on the particle's world line, so that the 
right-hand side of Eq. (p^q) is only formally defined; this expression must be regularized and the divergences must be 
removed. This renormalization of the self-force was carefully performed by Mino et al. and Quinn & Wald ]27| ] 
for the case of vacuum spacetimes. Fortunately, the extension of their work to nonvacuum spacetimcs is trivial, and 
it leads to an identical result. The reason is that the renormalization procedure is based on a local analysis that is 
sensitive only to the immediate vicinity of the world line. So long as the particle mov es in a region of spacetime that 
is empty of matter, which we assume here, the structure of the divergent terms in Eq. (|6.8| ) is the same whether o r no t 
matter is present somewhere in the spacetime. After renormalization, the self-force is found to be given by Eq. ( [1.3] ), 
but it is now expressed in terms of a Green's function that is sensitive to the presence of matter. 

We conclude that in a spacetime that contains matter, the gravitational self-force is given by its usual expression 



(1.3), but that it now involves a Green's function that incorporates information about the matter distribution; this 



information is encoded in the tensor S^avp defined by Eq. ( p. 14 ) 



We note that the renormalization procedure removes the 0(m) part of the self-acceleration, which now admits the 
expansion 

/« lf /m = O(mM) + 0(mM 2 ) + • • • . (6.9) 



in the background spacetime. Instead, 



Using this together with Eqs. (6.3) and (6.8), we find that the particle's stress-energy tensor is no longer conserved 
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r^.p = 0{m 2 M) + 0(m 2 M 2 ) + ■■■. (6.10) 

This is as it should be, because the motion is no longer geodesic in the background spacetime. In a formal sense, r"' 3 
is now conserved in a spacetime with metric gap ~t~ ^fap- 



C. Matter- mediated force 

The existence of additional terms on the right-hand side of Eq. ([T^) follows directly from the fact that the stress- 
energy tensor T a ^[g] describing the matter distribution depends on the metric. Because inserting the particle in the 
spacetime has the effect of shifting the metric from g a p to g a p + J a p, there must be a corresponding shift in the 
stress-energy tensor. Physically, this corresponds to the facts that (i) the particle raises a tide on the star and induces 
internal motions within the fluid, and (ii) the particle sets the star in a small motion around the system's center of 
mass. Both effects are incorporated in the shift in the stress-energy tensor; as we shall see, however, the bulk motion 
of the star is much more important for our purposes. We let 

ST 01 ? = T af3 [g + j] -T af3 [g] (6.11) 

be the perturbation of the star's stress-energy tensor created by the presence of the particle. Below we will use 
T"' 3 = T a/3 [g + 7] to designate the perturbed stress-energy tensor of the moving star, and T Q/3 = T a ^[g] to represent 
the background values — the star at rest. Because T af3 = 0{M) + 0(M 2 ) + ■ ■ ■ and j a p can be expanded as in 
Eq. (|6.7D, we have that 



ST a/3 = 0(mM) + 0{mM 2 ) + ■ ■ ■ . (6.12) 

The perturbation in the star's stress-energy tensor must appear, together with r"^, on the right-hand side of 
the Einstein field equations. It will contribute an additional metric perturbation 5g a p, over and above the original 
perturbation ^ a p directly associated with the particle. Physically, the new perturbation represents the difference 
between the gravitational field of the moving star and the background field of the star at rest. Because both r Q ^ and 
ST a ° are quantities of the first order in m, the (linear) metric perturbations can simply be added; the full metric is 
therefore g a/3 + Sg a p + j a p. 

Imposing the Lorentz-gauge condition separately on both 5g a p and j a p, we have that 



5g^(x) = 4 J G a l, v ,{x,x')bT^\x')^d A x'. (6.13) 
Using the expansions ( |6.6| ) and ( 6.12| ), we see that Eq. ( 6.1 3| ) implies 



6g af3 = 0(mM)+0{mM 2 ) + ---. (6.14) 

The matter- mediated force corresponds to 5g a p acting on the particle. In the absence of j a p, this action can be 
described by the statement that the particle must move on a geodesic of the metric g* a ^ = g a p + Sg a p. Physically, this 
means that the particle is now subjected to the gravitational field of the moving star. Mathematically, the equations of 
motion would take the appearance of Eq. ( |6.8[ ), but with 5g a p taking the place of 70-/3. There is no need to renormalize 
this equation: Sg a p is smooth on the world line. 

Our conclusion at this stage is that in a spacetime that contains matter, the equations of motion of a massive 



particle are given by Eq. (6.1), with a self-force given by Eq. (1.2) and a matter-mediated force given by 



f&* = ~ 7 (V„ ;/J + Sg% a - 5g a * + u^Sg a p,^)u^. (6.15) 



According to Eq. ( 6.14 ) 



f*Jm = 0(mM) + 0{mM 2 ) + ■■■, (6.16) 

which should be compared with Eq. ([3^). We see that in a weak-field approximation, both forces contribute terms 
of order 0(mM) to the particle's acceleration. 

We have successfully extended the Mino-Sasaki-Tanaka-Quinn-Wald equations of motion to spacetimes containing 
matter. We emphasize that this extension does not rely on any weak-field assumption regarding the background 



spacetime. In the following subsections we compute the matter-mediated force of Eq. (6.15); this discussion will rely 
on a weak-field approximation. 
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D. Determination of T? 



The matter-mediated force of Eq. ( 6.1 5| ) involves the metric perturbation S g a g w hich is computed from 8T a @ , the 
perturbation in the star's stress-energy tensor; the conversion is given by Eq. ( 6.13 ). We have already introduced the 
notation T" 13 — T Q/3 + ST a/3 for the stress-energy tensor of the moving star. 

We now need to gather information about 5T a @ . Our main source will be the conservation equations for the full 
stress-energy tensor T a/3 + 5T af3 + r 1 



a/3 



a/3 



a/3' 



g a /3, <5r Q a the corrections generated by 



in a spacetime with metric g a p + Sg a g + lap- 

a/3 

i a p, and C^ a p those generated by j a p. For example, 



This metric comes with a 



connection r^„ , a + ST^ a g + C^gj in which T^ al3 are the Christoffel symbols constructed from the background metric 



1 



C % = 9 (Vaj/9 + 7 '/3;a ~ W) • 



(6.17) 



By virtue of Eqs. ( |6.7|) and (|6.14| ), we have 



C* M Q/3 = 0{m) + 0(mM) + 0(mM 2 ) + ■■■ 



and 



(6.18) 



ar* = O(mM) + 0(mAf 2 ) + ■ 



(6.19) 



We also have T^ aj3 = O(M) + 0(M 2 ) H and we shall use the expansions of Eqs. (§T|) and ( |6.12| ). 

The conservation equations can be expressed as 

T a/3 + r a/3 + grpafi + 2 c^ T 0)fi + 2T {a p ST^ 1 + 2C {a J8T^ 1 + t 0) ") + 28Y {a JT^ + 5T^ + t^) 



(6.20) 



The first term vanishes because the background stress-energy tensor is conserved in the background spacetime. The 
second term can be ignored, because by virtue of Eq. ( 6.10| ), it is of second order in m. The third of fourth terms 
have the expansions 0(mM) + 0(mM 2 ) + ■ ■ ■. Among the remaining terms we have (in schematic notation) TST = 
0(mM 2 ) + ■■-, CST = 0{m 2 M) + ■ ■ ■, Ct = Q(m 2 ) + 0(m 2 M) +■■-, 8TT = 0(mM 2 ) +■■-, STST = 0(m 2 M 2 ) +■■■, 
and STt = 0(m M) + • • •. We simplify Eq. (3.20) by discarding all terms that are second order in m. Incorporating 
at this stage our weak-field approximation, we also neglect terms that are second order in M. The remaining terms 
are all 0(mM) and we obtain 



0(mM 2 ). 



(6.21) 



This equation ( partia lly) determines <5T Q/3 in terms of known quantities. Because we have only four equations for ten 
unknowns, Eq. ( 6.21 ) cannot be used without first adopting a specific phenomenology for the matter distribution. We 
shall come back to this point in Sec. VI F. 

It is important to notice that because Eq. ( 3.21 ) has been truncated to order mM, only the 0(m) part of 
is invol ved in the determinatio n of 8 T al3 

cf. Eq. (U). By virtue of Eqs. (|J), (|3~30|) , and (|6.5|),lhls is given by 



cf. Eq. (6.18). This can be obtained from the 0{m) part of 7 Q ^ 



7flft(^)=4 



F a/3 



(*■ 



X I' X ) j3„/ 



IX — X 



(6.22) 



where f Q/3 = r Q/3 — \{rj^T^ l ') r q al3 is the particle's "trace- reversed" stress-energy tensor. The subscript "flat" indicates 
that the metric perturbation is calculated as if the particle were moving in flat spacetime (on a Newtonian trajectory 
around the star); the right-hand side of E q. (|6.22 ) involves only the 0(m) part of r Q/3 — cf. Eq. (6.4). 
Another consequence of truncating Eq. ( 6.21| ) is that it can be re-expressed as 



-,a/3 



C a T^ + C%T^ = 0(mM 2 ), 



Pi 

nap 



(6.23) 



in terms of the star's perturbed stress-energy tensor: T"^ = T a @ + ST al3 . Here we use the facts that (i) T"^ is 



actually zero for the weak-field spacetimes considered in this paper — T .„ is nonzero but is a quantity of order M 
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— and (ii) the difference between CT„ and CT is of second order in m. Equation ( |3.25 ) implies that is conserved 
in a spacetime with metric g^ rtlclc = r\ a p + up to terms of order mM 2 . We will refer to this observation below. 

Let us n ow ta ke stock. In order to compute the matter-mediated force of Eq. (6.15) we must first calculate 5g a p 
using Eq. ( 6.13 ). We therefore need to find ST a/3 , the perturbation of th e stress-energy tensor associated with the 
star's motion around the center of mass. For thi s we shall use Eq. ( |3.21 ), after adopting a specific phenomenology 
for the matter distribution (Sec. VI F). But Eq. ( 6.21 ) involves also the connection coefficients C M Q ^ computed from 
lad - T° calculate this is our next task. 



E. Computation of 7^ 



The stress-energy tensor of the particle is given by Eq. (6/2). For the purpose of computing we let the particle 
move in flat spacetime, so that the metric used in Eq. (3.2) to define proper time r and the determinant g is r] a p, the 
metric of flat spacetime. After changing the variable of integration to z° and integrating over 5{t — zq), we obtain 



r a/3 (t, x) = midt/dr^v 13 5 3 (x - z), 



(6.24) 



in which z(t) represents the trajectory of the particle, v a — (1, v) with v(t) = dz/dt, and dt/dr = (1 — v 2 )^ 1 / 2 . We 
recall that although the calculation is carried out in flat spacetime, the particle is actually moving on a geodesic of 
gap; its motion is governed by Newton's equation, 



dv 



= g = -V$, 



(6.25) 



where the right-hand side is evaluated at x = z(t). 

It will prove sufficient for our purposes to evaluate T a ^ to second order in v, the particle's velocity. After trace 
reversal, we obtain 



-00 



-Oa 



m 
~2 



l + ^ 2 + 0(« 4 ) 



hix - z), 



mv a [l + 0(v 2 )] 5 3 {x - z) 
l-\v 2 + 0{v±) 



(6.26) 



_ {5ab 



+ 2v a v b [l + 0(v 2 )] \S 3 (x-z) 



To be consistent we also Taylor-expand the right-hand side of Eq. (3.22) about the current time t. This gives 



7 fl tt(^)=4 



T af3 (t, X') 

\x-x'\ 



d 3 x' - 4 



F a(3 



(t,x')d A x ' + 2 



dt 2 



f a(3 (t,x')\x-x'\d 3 x' 



It is now straightforward to substitute Eqs. (6.26) into Eq. (3.27) and perform the integrations. For example, 

2mA 



7flat 



2m m +m df 2Mx - zl 



(6.27) 



(6.28) 



where A = 1 + ^v 2 + 0(v 4 ). This result simplifies considerably: When the time derivative acts on v, it generates 
via Eq. ( 6. 25] ) a term of order g ~ M/\z\ 2 that can be neglected, because we are interested only in the 0(m) part of 
"fad- We are therefore left with evaluating d 2 \x — z\/dt 2 , which is a simple task. Subjecting the other components to 
similar manipulations, we arrive at 



where n = (x — z)/\x — z\. 



2m 
\x-z\ 
4mv a 
\x-z\ 

2m 
\x - z\ 



l + 2v 2 



1 



(n • v) 2 + 0(v 4 ) 



O(mM), 



[l + 0(i; 2 )] +0{mM), 



(6.29) 



§ ab 



l~-(n-v) 2 



+ 2v a v b + 0(v 4 ) \+0{mM), 
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F. Determination of T, (continued) 



As we have remarked at the end of Sec. VI D, the determination of T"^ via the conservation equation ( 6.23 ) is 
possible only if we adopt a specific phenomenology for the matter distribution. We shall make the simplest choice, and 
describe the star as a particle of mass M moving on a trajectory Z(t) with a velocity V(t) = dZ/dt. The pointlike 
nature of the moving star will not be a problem here: First, for the purposes of this calculation the star is not 
subjected to its own field because, as we have pointed out, it moves in a spacetime with metric <7^t rtlcle = r) a p + 7^; 
second, the singularity in 5g a p is far removed fr om t he particle's world line and is therefore not an obstacle in the 
calculation of the matter-mediated force of Eq. (6.15). Notice that in the following developments we shall formally 
treat both Z and V as quantities of the first order in to. 

Following the steps leading to Eq. ( |6.24 ) , we find that the stress-energy tensor of the moving star takes the form 



M- 



dt/dr 



V ^/particle 



V a V 6 3 (x - z), 



(6.30) 



where V a = {1,V) and dt/dr = (-gP^^WV)- 1 / 2 . After discarding terms of order to 2 , this reduces to T°° = 
MS 3 (x - Z) + M( 7 g° t - ±5 ab j^ t )6 3 (x), T° Q = MV a 6 3 {x), and Tf = 0(m 2 M). The perturbation 5T a $ can be 
obtained from these results by removing the background values, T 00 = M5 3 {x) and T 0a = T ab = 0. 

Substitution of our expressions for 5T af3 into Eq. ( 6.21 ) leads to equations of motion for the star. The spatial 
components of the conservation equation imply 



dV a 1 d 



At 



2<9r 



d_ 

dt 



j° R : t + o(tom), 



(6.31) 



in which the right-hand side is evaluated at x = 0. We find that the time component of the conservation equation 
holds as an identity, but this confirms the necessity of including y/— <7 pa rticic in our expression for the stress-energy 
tensor. 



It is a simple task to evalu ate th e partial derivatives of Eq. ( 6.31 ) starting from the expressions listed in Eq. ( 3.29 ) 
in this procedure we use Eq. ( 6. 25] ) to dismiss dv/dt as a term of order M/\z\ 2 . After evaluating the results at x = 
we obtain 



dV 
~~dt 



1 + 2v 2 - -(£ • vf + 0(v 4 



3(2 • v)v } + O(mM), 



(6.32) 



where z = z/\z\. The mot ion of the star is now determined to a degree of accuracy sufficient for our remaining 
calculations. Equation (|6.32| ) confirms that V, and therefore Z, are to be treated as quantities of order to. We note 
that in th e present context in which the star has become a point mass, the motion of the particle is still governed by 
Eq. ( 6.25 ), but that the gravitational field g now takes the explicit form 



71/ 



(6.33) 



Equation ( 6.25 ) determines the vectors z{t) and v(t) = dz/dt that appear in Eq. ( 6.32 ). 

The perturbed stress-energy tensor of the moving star is now completely determined. Its final expression is obtained 
by substituting the potentials of Eq. (3.29) into the results given before. This yields 



T™(t,x)=M8 3 (x-Z) 



T» a (t,x) = MV a S 3 (x), 



raM 



1 - 2v 2 - hz ■ v) 2 + 0(v 4 ) 



S 3 (x) + 0(mM 2 ), 



(6.34) 



and T. 



ah 



0(m 2 M). 



G. Computation of 



a/3 



The metric g*^ = g a /3 + 5g a p was introduced in Sec. VI C, in the paragraph following Eq. (6.14). It represents 



the gravitational field of the moving star, and the action of the matter- mediated force of Eq. ( |6.15 ) can be described 
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by the statement that the particle must move on a geodesic of this metric. In the next subsection we will use this 
statement to produce the equations of motion of the particle under the sole action of the matter-mediated force. In 
this subsection we calculate g* a a. 

T he me tric of the moving star differs from g a p, the metric of the star at rest, by Sg a p, which can be computed from 
Eq. (3.13). Keeping only the O(mM ) part of this equation, which involves the O(l) part of the gravitational Green's 
function — cf. Eqs. (3J3) and ( 3.3Cl|) , we write this as 



5g a *{t,x) 



5T al3 (t- \x-x'\,x') o . , 9x 
jr— - — - d 3 x' + 0(mM 2 ), 



(6.35) 



where ST af3 = 6T af3 — ^{r} tlv 5T fiv )rf x ^ . At the same time, we recall from Sec. II A that g a p differs from r) a p by 



quantity h a £ 



-2$Xq/3 that can expressed as 



h af3 (t,x) 



T al3 (t 



\x — x " , as 



d 3 x' 



(6.36) 



where T a P is the trace-reversed stress-energy tensor of the background spacetime. Adding Eqs. ( 6.35| ) and ( 6.36| ), we 
find that 



where h* af3 = h a p + bg afi is given by 



9 *a(i = Va0 + h* af3 , 



ff(t-\x-x%x') 



as — x 



d 6 x' + 0{mM 2 ), 



(6.37) 



(6.38) 



and where = T" /3 - \{^ V T^ )^ fi can be calculated from Eq. ( |Jg) . To first order in m, this gives T°° = |T°°, 



f° a = T° a , and Tf = \5 ah T2°. 



We shall describe in some detail the computation of We first Taylor-expand T°°(i — |as — x'\,x') about the 
current time t and express h® as 



h°° = 2 



~lw J T *°^ x ')\ x - X '\ 2 d 3 x' + --- + 0(mM 2 ). 



(6.39) 



We will see that truncating the series after four terms gives us sufficient accu racy: it is very important, however, to 
keep the last term involving three time derivatives. After substitution of Eq. (3.34), we obtain 



,00 



2M 
\x-Z\ 



Q2 i g3 

M ^- Z '-3 M ^- Z ' 2 



2mM A n ,J A , x . . d 2 A 1 , ,. .„ d 3 A 
— — — — + 2mM— — - mM as — — + -mM \x\ 2 — — 
x z dt z 1 l dt 2 z 3 11 dt 3 z 



0(mM 2 ), 



(6.40) 



where A = 1 - 2v 2 - \{z ■ v) 2 + OM ), It is a straightforward task to evaluate the time derivatives contained in 
Eq. (3.40); for this we involve Eq. ( |3.32 ) and we make su re to discard all terms that are not first-order in both m and 
M. For the first line (fl) on the right-hand side of Eq. ( 6.4C ) we obtain 



2M x ■ z r , 9n1 2 do r . 9n1 . 9n 

fl = r — - mM—r—p. 1 + O - -mil ~ 1 + 0{v 2 )] + 0(mM 2 ), 

\x — Z\ as z r 3 dt 



(6.41) 



where g is defined by Eq . (6.33). In Eq. ( 3.40] ), the first set of "• • •" terms represent a correction of order w 4 to the 
second term of Eq . ( 3.41 ); these can be safely ignored. The last member of Eq. ( 6.41 ) comes from the third-derivative 
term in Eq. ( 6.40| ); this will be seen to give rise to a radiation-damping term in the matter-mediated force. For the 
second line (sl) on the right-hand side of Eq. ( 6.40 ) we find 



2 777 Z ■ 1) 

SL = -T-rr-r [l + 0{v 2 )} - 2mAf— [l + 0{v 2 )] + 0{mM 2 ). 
\x\\z\ \z\ 



(6.42) 
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Here the third-derivative term in Eq. ( |6.40 ) is seen to give rise to the 0{v 2 ) correction in the second member of 
Eq. ( 6.42 ), and the "• • •" terms are responsible for 0(v A ) corrections in the first member. 
Gathering the results, we arrive at 



hf{t,x) = 



2M 
\x-Z\ 



mM 



x z 



X ■ z 



2 

— —mx 



dq „ , r z 
-f- - 2mA/— 
dt \z 



(9(2pn,toM 2 ), 



(6.43) 



in which 0(2pn) regro ups a ll the 0(v 2 ) correction terms displayed in Eqs. ( |6.41 ) and ( |6.42 ). Wc note that on the 
right-hand side of Eq. (6.43), the first member is of order <!> and is therefore designated as "Newtonian" . The second 
member is of order m/\z\ relative to <&, and it represents a post-Newtonian (Ipn) correction |^] . Because they involve 
an additional power of v, the last two members are 1.5pn corrections. 

The calculation of h® a involves similar steps, but the computations are substantially simpler. We simply quote the 
result: 



hl a {t,x) 



4MV a 



x 



+ Amg a + 0(1. 5pn, mM 2 ). 



(6.44) 



Relative to the leading, "Newtonian" term in the first term on the right-hand side of Eq. ( 3.44 ) represents a 
0.5pn correction, while the second term is a Ipn correction. Because /i° a couples to the particle's velocity v a , these 
post-Newtonian labels should actually be promoted to Ipn and 1.5PN, respectively; the error term is then 2PN and 
h° a has been computed to a sufficient degree of accuracy. 



The metric of the moving star is now completely determined. We have, from Eq. (6.37) 



Soo 



,oo\ 



0{mM 2 



90a 



hl a + 0(mM 2 ) 7 g* ab = S ab (l + hf) + 0(mM 2 ), 



(6.45) 



with h®° and h® a given by Eqs. ( |6 . 43 ) and ( |3.44 ), respectively. In those expressions, z(t) represents the Newtonian 
motion of the particle around the fixed star, v(t) = dz/dt, and dv/dt = g = —Mz/\z\ 3 . On the other hand, Z(t) 
represents the orde r-m m otion of the star in the gravitational field of the particle, V(t) = dZ/dt, and dV /dt is 
determined by Eq. ( |6.32 ). 



H. Motion of the particle under the matter-mediated force 



In this subsection we calculate the motion of the particle under the sole action of the matter-mediated force; 
we neglect for now the action of the self-force, which will be incorporated later. In this subsection, ther efore , we 
shall write the equations of motion as u°„/ = f^ m - The force could be calculated directly from Eq. ( |6.15| ) by 
substituting 8g a p obtained from Eq. ( 6.45| ). It is easier, however, to proceed in the following way. First, we recognize 
that the equations of motion are equivalent to the statement that the particle moves on a geodesic of the metric 
9af3 = 5 'a/3 + ^9a(3- Second, we generate the equations of motion by constructing a suitable Lagrangian function 
L(x, x) which we then substitute into the Euler-Lagrange equations. The practical advantage of this method is 
that it allows us to straightforwardly switch from a world line parameterized by r, the particle's proper time on the 
background spacetime, to one parameterized by t, the time coordinate. In the following we will denote the updated 
trajectory of the particle by x{t) and its updated velocity by x(t) — dx/dt. The Newtonian values will still be denoted 
z(t) and v (t). Because the matter-mediated force provides a correction of order m to the Newtonian motion, we have 
x(t) = z(t) + 0(m) and x(t) = v(t) + 0(m). In the next subsection we will abolish this distinction of notation and 
use z and v to refer to the updated trajectory. 

We take the Lagrangian function to be 



L(x, x) — 1 



-9W 



(6.46) 



where x a = (1, is); as we have indicated, we use t as a parameter on the world line. We substitute the metric of 



Eq. (6.45) into Eq. (3.4£), and to simplify we expand the square root to second post-Newtonian (2pn) order; in this 



procedure we also discard terms that are second order in M. The result is 



L 



1 



1 



-h 



oo 



1 



h° a i a 



h°°x 2 + 0(2pn,M 2 ). 



(6.47) 



After substitution of Eqs. (6.43) and (6.44), we notice that the Lagrangian contains a term —mMz ■ v/\z\ 3 . This is 
a function of time only which does not affect the equations of motion; we therefore remove it from the Lagrangian. 



30 



The Lagrangian also contains a term —Amx ■ g. This differs from 4mx ■ dg/dt by a total derivative —Amd(x ■ g)/dt 
that can also be deleted from the Lagrangian. The final result is an effective Lagrangian 



1 



M 



1 



mM 



1 



2d 



AM 



V 



11 



dg_ 

dt 



0(2pn,M 2 ) (6.48) 



that is ready to be substituted into the Euler-Lagrange equations. 
The first member of the Euler-Lagrange equation is 



d dL 



cir 



dt dx a 



A 



ab 



dx b 
~dt 



3M 



(x- 



3M 
UI3 ( x 



V + x- Z)x a 



AmM 



(6.49) 



where A ab = (1 + \x 2 + 3M/\x - Z\)S ab + x a x . To obtain this we have used Eq. ( |6.32| ), appropriately truncated at 
Newtonian order, to evaluate dV a /dt. The second member is 



dL. 



cir 



dx a 



-Mi 



. x a - Z a 



Z\* 



1 



mM 



1 



x ■ z 



mM 
2\x\\z\ 



AM r. t t\ o, 11 d 9 a 
r-Tg (x ■ V)x a + —m^- 



x 



dt 



(6.50) 



The equations of motion are obtained by equating Eq. ( 3.49 ) to Eq. ( |6.50 ) and solving for dx a /dt. This requires the 
inversion of A ab , which can be performed approximately to first post-Newtonian order. To simplify our expressions 
we replace z a by x a + 0(m) in the terms proportional to mM, and we also replace x by x — Z = x + 0(m) in terms 
that are already linear in to. (It is of course important that we make these substitutions after varying the Lagrangian, 
and not before.) 
The final result is 



dx a 
~dT 



M 



Z\* 



l + x 1 



5m 



Ax ■ V 



Z a ) ~{x-Z)- {Ax - 3V) x a + A(x-Z)-x V a 



li m ^! +0 (2 PNj M 2 ). 
3 dt v ; 



(6.51) 



We recall that Z(t) = 0{m) represe nts t he trajectory of the moving star, that V(t) = dZ/dt is determin ed by 
Eq. ( 6.32 ), and that g, given by Eq. ( |6.33j ), is the Newtonian gravitational field of the star at rest. Equation ( 6.5l| ) 



governs the motion of the particle under the sole action of the matter-mediated force. We see that the last term, 
which represents a 1.5pn correction to the Newtonian equations of motion, gives rise to radiation damping. This 
term, however, is precisely canceled out by the gravitational self-force of Eq. ( |5.3l| ). We conclude that the motion of 
the particle, under the combined action of the self-force and the matter-mediated force, is conservative at this level 
of approximation. 

I. Summary and comparison with post-Newtonian theory 

This is the end of the line, and we better summarize. A particle of mass to moves in the gravitational field of a 
star of mass M. The particle's trajectory is denoted z{t) and its velocity is v(t) = dz/dt. The star is also moving, 
on a trajectory Z(t) with a velocity V(t) = dZ /dt. We let p = z — Z be the separation between the two objects, 
and we use the notation p = \z — Z\. We let g = —Mz/\z\ 3 be the Newtonian gravitational field of the star at rest, 
g* = —Mp/p 3 is the gravi tatio nal field of the moving star, and 5g = g* — g is the difference between the two. The 
spatial components of Eq. (6.1) then give 



./self fn 



dv 

dt a to m 

for the particle's equations of motion. The gravitational self-force was computed in Sec. V C and is given by 

/ S cif 11 dg 



TO- , 

3 dt 



<3(2pn, ikT). 



The matter-mediated force was computed in the preceding subsection and is given by 



Sg 



5to 
P 



Av ■ V 



M 

"^3 



p- (Av - 3V)v - A{p ■ v)V 



11 dg 
— to — 

3 dt 



0(2pn,M 2 ). 



(6.52) 



(6.53) 



(6.54) 
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The final form of the equations of motion is obtained by combining Eqs. ( |6.52 )-( 5.54). The result is 



dv 
~dt 



M 
7 



, 5m 
1 + v 2 



4v ■ V )p - p ■ (4v - 3V)v + 4(p ■ v)V 



C>(2pn,M 2 ), 



(6.55) 



and this contains no trace of a radiation-reaction force. The motion of the star, on the other hand, is governed by 
Eq. (6.32), which we rewrite as 



dV 

dt 



in 
1i 



P 



1 + 2v 2 



3 ( p ■ v 

P 



p - 3(p • v)v } + 0(2pn, mM). 



(6.56) 



In going from Eq. (6.32) to Eq. (3.56) we have allowed for the fact that p=z — Z = z + 0{m). 

The reader familiar with post-Newtonian theory (as presented, for example, in Ref. [|60|) will have undoubtedly 
noticed a similarity between the calculations presented here and the standard post-Newtonian treatment of a system 
of point masses. While similarities are surely present, the differences are important. The most noticeable difference 
is concerned with self-field effects, which are sometimes swept under the rug in the post-Newtonian treatment — the 
mass of each body is simply renormalized every time an infinite self-field appears in the equations. By contrast, the 
computations presented in this section are completely free of infinities — those have been taken care of, once and for 
all, during th e con structio n of the self- force. 

Equations (6.55) and (3.56) should be compared with the equations derived from a standard post-Newtonian 
treatment — the Einstein- Infcld-Hoffmann equations of motion |j ] . These are 54 



dv 
~dt 



AM 
P 



5m 
P 



+ v 2 + 2V 2 



2\ p 



p - p- (4w - 3V)(v - V) } + 0(2pn) 



.57) 



and 




1 - 



4m 



5M 
P 



V 1 



2v 2 - 4v • V - 



3fp- 



p + p- (4V-3«)(«- V) } +0(2pn). (6.58) 



If we eliminate all terms of order m 2 and M 2 from Eq. (6.57), we recover Eq. Q6.55 ). If, on the other hand, we remove 



all terms of order m 2 and mM from Eq. ( 6.58 ), we recover Eq. ( 6.56 ). Our results are therefore perfectly consistent 



with the standard post-Newtonian treatment. 
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APPENDIX A: HADAMARD FORM OF THE SCALAR GREEN'S FUNCTION 



In this appendix we continue our discussion of the two-point function A{x,x'). While the results derived here are 
not required for the computation of the self-forces, they are interesting in their own right. Our main objective is to 
show that the scalar Green's function computed in Sec. Ill B can be cast in the Hadamard form |29 



G{x, x') = 9(x, x') u(x, x')S(a) + v(x, x')6(-cr) 



(Al) 



Here, a(x, x') is Synge's world function |n]], equal to one-half the squared geodesic distance between x and x'; u(x, x') 
and v{x,x') are two-point functions that are smooth at a = 0, and 0(x,x') is a time-ordering function, equal to 1 
if x is in the causal future of x' , and zero otherwise. The calculations presented in this Appendix rely on methods 
introduced by Thorne and Kovacs [|62|. 
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1. A(x,x') on the light cone 



We go back to Eq. (4.6) and seek to evaluate A(x 1 x') in the situation where the points x and x' are linked by a 
null geodesic of the flat-spacetime background. In this situation At = R, or s = e. (Re call that At = t — t' and 
R = \x — x'\.) This means that r\ l and rj 2 are both zero, and that the vector of Eq. ( |4.5| ) reduces to r\ — ^cos9R. 
(Recall that the third axis is oriented along R = x — x' .) Integrating over <fi, we find that Eq. (4.6) becomes 

A(x,x') = I &(xq + rj) dcosd. After changing the variable of integration to A = i(cos# + 1), this becomes 



Ai c (x,x') 



where 



*(€(A))dA, 



(A2) 



£(A) = x' + A(sb - x') 



(A3) 



is a vector that interpolates linearly between x' (when A = 0) and x (when A = 1). In Eq. (AS) we have indicated 
with the label "lc" that A(x, x') is evaluated on the light cone. We see that A\ c (x, x') is the average of the Newtonian 
potential $ over the straight line joining x and x' . 

Up to now, our view of the two-point function A\ c (x,x') has been that it is the restriction of A(x,x') on the light 
cone: the points x and x 1 must be related by a null geodesic of the background Minkowski spacetime. The right-hand 
side of Eq. ( |A2| ), however, stays meaningful even for points that are not so related, and we take this opportunity to 
extend the definition of A\ c (x,x') to arbitrarily related points. This will be our new point of view: The two-point 
function A\ c {x,x') shall be defined for arbitrary points x and x' by Eq. (A2). In the case where x and x' are null 
related, then A\ c {x, x') — A(x,x'). 



We now show that A\ c (x, x'), as defined for arbitrary points by Eq. (A2), is intimately related to the world function 
a(x,x'), defined by (62] 



a(x, x ) 



1 tt*^ *\ 

9a ^ ] ^x^x dA ' 



(A4) 



where £ Q (A) describes the (unique) geodesic linking the points x and x'; A is an affine parameter scaled in such a way 
that £"(0) = x a and £ Q (1) = x a . The world function is to be calculated with the metric of Eq. (|2.3|). 



In Eq. (A.4), the geodesic 7 differs from a straight line by a quantity of order $. Because Eq. (A4) is an action 
principle for the geodesic equation, this error of the first order in the specification of the curve produces an error of 
the second order in a. To calculate a to first order in $, it is sufficient to take 7 to be a straight line . The appropriate 
relations are then £ a (A) = x a + X(x a — x a ). Substituting also g a/ 3 = r) a p — 2$(£)x a /3 in Eq. (A4), we obtain 



<j(x, x 1 ) = a 6at (x, x') - (At 2 + R 2 ) f 1 $(C(A)) dX, 

Jo 



(A5) 



where cr flat (a:, x') = ii_-At 2 + R 2 ). 

Comparing Eqs. (A2) and (A5), we arrive at 



a(x, x') = CT flat (x, x') - (At 2 + R 2 ) A lc (x, x'). 



(A6) 



We see that the two-point function A\ c {x, x') determines by how much the squared geodesic interval a differs from its 
flat-spacetime value. 



2. van Vleck determinant 

The scalarized van Vleck determinant | ^l| , |63t plays a fundamental role in the theory of Green's functions in curved 
spacetime. This is defined by 

v /g(x)g{x') 
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where a(x,x') is the world function defined in Eq. (A4), and g(x) is the metric determinant evaluated at x. The 
general theory predicts that in Eq. ( Al), u(x, x') is the square root of the van Vleck determinant. We will verify this 
result with our we ak-c urvature scalar Green's function. In order to do so we must compute A(x, x'). 
We rewrite Eq. (A£) in the form a = cr flat — SI, where 



n{x,x') = X ^(x-xY{x-x'yA lc {x,x'). (A8) 



A straightforward computation, using y / —g(x) = 1 — 2$(a;) and d a /3> cr flat = — r] a p, reveals that to first order in the 
Newtonian potential, the van Vleck determinant is given by 

A = 1 + 2$(x) + 2$(x') + rf^d^-n. (A9) 

The last term is evaluated from Eq. ( |A8|) , and we obtain 

V^d^tt = (At 2 + R 2 )5 ab d ab ,A lc - 2(x - x') a (d a - d a ,)A lc - AA ic , (A10) 

where we have used the fact that A\ c (x,x'), as defined by Eqs. ( |A^ ) and (|A3|), depends only on the spatial variables 
x and x' . Calculating the derivatives is easy, and we get d a A\ c — J A$, a dX, d a >A\ c — J(l — A)<I>, a dX, and d a b'A\ c = 
J A(l — A)$ !a fc dX. Inside the integrals, the derivatives of $ are taken with respect to £, and the vario us fac tors 



involving A come from differentiating this vector with respect to either x or x' . Substituting this into Eq. ( A10 ), we 
have 

T] al3 d al3l n=(At 2 +R 2 ) f X(l-X)V 2 <S>dX-2(x-x') a f (2A - 1)$ „ dX - 4 / $dA. (All) 
Jo Jo Jo 

In the first integral we replace V 2 $ by Anp. In the second integral we replace & y a(x — x') a by d^/dX, which allows 
us to integrate by parts. After simplification, we arrive at 

T] al3 d a f3 / n = 4ir(At 2 + R 2 ) [ A(l - X)pdX - 2$(x) - 2<$>(x'). (A12) 

Jo 

Substituting this into Eq. ( |A9| ), we obtain our final expression for the van Vleck determinant: 

A{x,x') = l + 4n(At 2 + R 2 ) f A(l - A)p(^(A)) dX. (A13) 

Jo 

We see that A(x,x') differs from 1 if and only if the straight line connecting the points x and x' passes through the 
matter distribution. 

3. A(x, x') near the light cone 

We have previously evaluated the two-point function A(x,x') in the case where x and x' are connected by a null 
geodesic of the background Minkowski spacetime. We now improve on this result, by allowing the parameter 

e = Vs 2 - e 2 = - yjAt 2 - R 2 (A14) 

to be nonzero. We will, however, assume that e<1, and evaluate Eq. ([O]) in this limit. This is the near- light-cone 
approximation. 



To second order in e, the vector r) of Eq. (4.5) reads r/ 1 = es'm6 cos0, rj 2 — £sin6*sin0, and rj 3 = e[l + \ieje) 2 + 
0(e 4 )] cos 9. We substitute this into &(xo + T]), expand to second order in e, and then integrate over the angles. After 
re- introducing the parameter A = ^(cos# + 1), we obtain 

A = A lc + 4ire 2 f A(l - X)pdX + ^- [ (2A - 1)<!>> Q dX - e 2 f A(l - A)$, a6 n a n b dX + 0(e 4 ). (A15) 
Jo 2e J Jq 

Here, A\ c is the restriction of A(x, x') on the light cone, the quantity given by Eq. ( |A2j ), and n is a unit vector po inting 
in the direction of x — x' . The quantities inside the integrals are evaluated at the point £(A) defined by Eq. ( AS ), and 
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the derivatives of $ are taken with respect to this vector. To get the first integral we have used Poisson's equation, 
V 2 $ = 47175. 

In the second and third integrals of Eq. (A1E), the derivatives of $ in the direction of h can be express ed a s 
derivatives with respect to A: d^/dX = 2e$ ja n a and d 2 $/dX 2 = (2e) 2 $ t abh a h b . After substitution into Eq. (A15) 
and an integration by parts on the last term, we find that these integrals cancel out. Our final expression for the 
two-point function is therefore 



A(x, x') = A lc (x, x') + Aire 2 / A(l - A)p(£(A)) dX + Q(e 4 ). 



(A16) 



We recognize here the same integral over the mass dens ity t hat appears in our previous expression for the van Vleck 
determinant, Eq. ( A1S ). This allows us to rewrite Eq. ( A16) as 



A(x, x') = A lc (x, x 1 ) + gR2 n [A(x, x') - 1] + 0(e 4 



(A17) 



In this calculation we have used the fact that the factor At 2 + R 2 appearing on the right-hand side of Eq. ( A13 ) is 
equal to 2R 2 [1 + 0(e 2 )}. 

Equations ( |A16| ) and ( A17 ) indicate that unless the straight line connecting x to x' intersects the matter distribu- 
tion, Ai c (x,x') makes an excellent approximation to A(x,x r ). 



4. Scalar Green's function 



We now have the necessary tools to cast the s calar Green 's fu nction calculated in Sec. Ill B into the Hadamard 
form (2^] displayed in Eq. (Al). Recall from Eq. ( [3.11 ) and ( |3.14 ) that our expression for the Green's function was 

G(x,x) = G aat (x,x r ) - 2d tt >A(x,x / ) ~ 2£,B(x,x'), (A18) 

where G flat (x,x') is the retarded Green's function of flat spacetime, given by Eq. ( |3.3| ). This can be re-expressed as 

G &at {x,x') = 0{x,x')6(<j aat ), (A19) 

where 9(x,x') is the time ordering function introduced in Eq. flAl[ ), and a &at (x,x') = |(— At 2 + R 2 ). 

To calculate the Green's function we need to take into account the fact that both A(x,x') and B(x,x') are zero if 
x and x' are spacelike related (in the flat-spacetime background), and if x' lies to the future of x. We express this as 



A(x,x') =9{x,x')9(~a fiat )A{x,x'), B(x, x') = 8(x, x')6(-a Rat ) B(x, x'), 



(A20) 



in which the 0-functions explicitly enforce the vanishing of the two-point functions except when x' is in the causal past 
of x. In what follows we will keep the time-ordering function implicit, and re-insert it at the end of the calculation. 
Differentiation of Eq. (A20) gives 



dtt>A = e(-a aat )d tt ,A - J(<T flat )(l + 2Atd t )A + 6'(a Rat ) At 2 A, 



(A21) 



where the prime on the <5-function indi cates differentiation with respect to cr flat . To calculate the time derivative of 
A in the second term, we use Eq. ( A16| ) which we copy as A — A\ c — a &at C + 0(cr flat2 ), where C — 2ir j A(l — X)pdX. 
Dropping all terms proportional to cr fiat , we find that the coefficient of the <5-function becomes — A\ c — xC, where 
X = At 2 + R 2 . Working now on the third term, we find that the coefficient of the differentiated ^-function is 
Ix^ic — (Aic + ^xC') fjfiat + 0(er flat2 ). The term that is linear in er flat can be transferred to the coefficient of the 
(5-function with the help of the distributional identity <t5'(<j) — —5{<j). The end result is 



d w A = 6(-<r &at )d w A- \ X C5{a iiat ) + \ X Ai c 5' (a flat ) . 



We substitute this into Eq. QA18p . With the help of Eqs. ( |A19| ) and ( |A20| ), we obtain 
G = (1 + X C)5{a VLat ) ~ xAcS'ia^) + 6»(~CT flat ) \-2d H >A - 2£B 



(A22) 



(A23) 
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We can re-express the first three terms in the form (1 + xC)5(cr flat — xAic)- From Eq. ( |A6| ) we recognize the new 
argument of the (^-function as the world function cr(x, %') of the weakly curved spa cetim e. And we recognize the factor 
1 + xC in front as the square root o f the van Vleck determinant, calculated in Eq. ( A13). Finally, because the quantity 
within the square brackets in Eq. (A22) is already of the first order in the Newtonian potential, we can safely replace 



a with a as the argument of the 0-function. The only remaining task is to re-insert the time-ordering function that 
has been left out of our expressions. 

Our conclusion is that the scalar Green's function of Eq. ( A1S| ) can indeed be cast in the form of Eq. (Al), with 
u — A 1 / 2 . For the weakly curved spacetimes considered in this paper, the world function cr(x, x') is worked out in 
Eq. (A6), the van Vleck determinant A(x,x') in Eq. (A13), and 



v(x, x') = -2d tt 'A(x, x') - 2£,B(x, x') 



(A24) 



is the tail part of the Green's function. Here, A(x, x') and B(x, x') are the two-point functions introduced in Sec. Ill 
A. 



APPENDIX B: TWO-POINT FUNCTIONS FOR THE POINT-MASS POTENTIAL 



In this Appendix we evaluate the two-point functions A(x, x 1 ) and B(x, x') for the special case $(cc) = —M/\x\. The 
results derived here were first obtained by DeWitt and DeWitt 0] ; we include this discussion here for completeness. 
We begin with Eq. (|4.6|), which we rewrite as 



A{x,xf) = ~ [ | — - — rdfl, 



(Bl) 



where rj(s, 8, </>) is the vector of Eq. (4.5) with s = ^At = h(t — t'), and »7o = — Xq = — ^{x + x'). The point rjo, at 
which the central mass is located (in a coordinate system centered at r\ — 0), will be represented by the ellipsoidal 
coordinates (sq, 9o, 4>o)- We recall that e = \R = h\x — x'\ is the cllipticity of the coordinate system. 
To evaluate the integral of Eq. (Bl) we invoke the addition theorem in ellipsoidal coordinates f&ifl , 



1 



-in 



\V-Vo\ 



Z=0 m=— I 



(i -my. 

(l + m)\ 



P l m (s</e)Qr(s > /e)Y l * m (d 0) c/>o)Y lm (0, </>), 



(B2) 



where s< (s>) is the lesser (greater) of s and so, and P™ and Q™ are associated Legendre polynomials [^5|. Substi- 
tuting this into Eq. ( |B1[ ), we find that the integration over the spherical harmonics Yi m (0,(j)) is zero unless / and m 
are both zero, and we obtain A = —(M/e)Po(s < /e)Qo(s > /e), or 



M 2 S> +R 



(B3) 



using the known forms for the Legendre functions of zeroth order. 

The hard part of the calculation resides with the computation of sq. Recalling that x and x' differ by a vector R 
pointing in the z direction, we write x = (x, y, z) 1 x' — (x, y, z — R), and we have r/o — (—x, —y, —z + ^R). This and 
Eq. (4.5) give us the equations 



sin f o cos (po = —x, yso 2 — e z sin 0$ sin <po — —y, sq cos 0q = —z + e, 
which must be solved for sq. It is not hard to show that this amounts to solving the quadratic 



so 



I2\ 2 i / 2 

T > S ° + 16^ 



0. 



(B4) 



(B5) 



where r = |x| and r' — \x'\. This finally gives so = |(r + r'). 

Going back to Eq. (B3), we see that we must distinguish between the cases At > r + r' , for which 2s> = At, and 
At < r + r' , for which 2s> = r + r' . This gives 



M 

A(x,x ) = — -0{At - R) 
R 



a ,\ -i r + r' + R At + R 

At) In + 9(At -r-r') In 



R 



At-R 



(B6) 



3G 



, and At = t — 
(O) and @ 



t'. For completeness we have re- inserted the step 
We see that the two-point function undergoes a 



where, we recall, r = \x\, r — \x \, R = \x — x 
function 6 (At — R) that was left implicit in Eqs 
change of behavior when At = r + r' . This time delay corresponds to a signal propagating with the speed of light 
from x to the center (at which the Newtonian potential is singular) and then on to x' . For shorter delays, A(x,x') 
is time-independent; for longer delays, A(x,x') depends explicitly on At. This sudden change of behavior makes the 
two-point function slightly suspicious: Although our calculations are based on the assumption that the spacetime is 
weakly curved everywhere, the change of behavior is dictated by a region of spacetime — the center — in which the 
Newtonian potential is decidedly not small. 

The computation of B(x,x') is quite simple for the potential $ = —M/r, for which the mass density is p(x) = 
M5(x). Substituting this into Eq. (4.2) gives 



M 



B{x,x') = — S(At-r-r'). 



(B7) 



Again we notice a suspicious dependence on the conditions at the center. 

Despite the notes of caution, the results der ived here for the spe cial c ase <£> = —M/r are in complete agreement with 
our general results of Sec. IV. First, Eq. (B7) is identical to Eq. ( 4.18 ), which was obtained as the leading term in a 
multipole expansion for B(x,x'). This indicates t hat in fact, our result for this two-point function is quite insensitive 
to the conditions near the center. Second, Eq. (B6) implies that near coincidence (R small and smaller than At), 
A(x, x 1 ) can be approximated by 



A(x, x) 



2M 
r + r' 



R 



r + r' 



0(R 4 ) 



(B8) 



which could b e rec ast in the form of Eq. (4.11). In particular, differentiation of Eq. (BS) with respect to x and x' 
confirms Eq. (4.13). In this case we have = (M/r 3 )(5 a b — 3r jt8 r b), with r a = x a /r. Third, Eq. ( pif ) implies that 
for long delays (At large and larger than R), we have the approximation 



At / \ 2M 

A(x,x') = - — 



1 1 / R 



O Ar 



(B9) 



which is evidently compatible with Eq. (4.14). This comparison between our general results of Sec. IV and those of 
this Appendix shows that any aspect of the two-point functions that might be sensitive to the strong-field portion of 
the spacetime near the center will not be involved in the computation of the self-forces. This statement provides a 
further validation of the work of DeWitt and DeWitt Em], which was entirely based on the special case $ = —M/r. 
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